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SUMMARY 


A  three-year  investigation  of  the  power  of  the  Euler  equations  in  the 
prediction  of  separated  flows  has  been  conducted.  The  preliminary  work  for 
this  study  was  conducted  from  6/84  to  6/85  under  AFOSR  Contract  F49620-84-C- 
0056  and  is  reported  upon  in  Part  I  of  this  report.  Here,  only  the  work 
conducted  in  the  last  two  years  of  this  research  project  will  be  discussed, 
except  when  background  material  from  Part  I  is  essential  to  the  completeness 
of  this  report. 

The  Euler  equations  are  solved  numerically  for  the  highly  vortical 
supersonic  flow  about  simple  geometries.  Two  sources  of  vorticity  are 
investigated;  the  first  is  the  flow  field  shock  system  and  the  second  is  the 
vorticity  shed  from  a  separating  boundary  layer.  In  the  case  of  shed 
vorticity,  the  surface  point  from  which  the  vorticity  is  shed  (i.e., 
separation  point)  is  determined  empirically.  Solutions  obtained  with  both 
sources  of  vorticity  are  studied  in  detail  and  the  relationship  between  the 
two  is  discussed.  The  complex  interaction  between  shock-vorticity-produced 
vortices  and  the  shock  system  itself  is  shown  to  result  in  interesting 
shock/vortex  patterns.  The  vorticity  shedding  model  developed  in  this  effort 
was  successfully  applied  to  flat-plate-delta-wing  and  simple  three-dimensional 
flows.  Less  success  was  achieved  in  applying  the  model  to  subsonic  flows.  A 
very  exciting  result  was  found  in  an  investigation  of  very  high  angle-of- 
attack  flows.  It  was  found  that  the  symmetric  solution  for  the  separated  flow 
about  cones  became  unstable  for  angles  of  attack  higher  than  two  and  a  half 
times  the  cone  half  angle.  At  these  angles  of  attack,  only  asymmetric 
solutions  exist,  as  noted  experimentally. 
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1.  INTRODUCTION 

Separation  and  the  formation  of  vortices  in  the  three-dimensional  flow 
about  wings  and  forebodies  at  high  angles  of  attack  constitute  an  important 
aspect  of  advanced  vehicle  design.  The  phenomenon  will  be  especially 
important  at  maneuver  conditions  for  the  next  generation  fighter  aircraft, 
since  the  flight  envelopes  of  these  vehicles  may  extend  to  high  angles  of 
attack  at  supersonic  speeds.  An  understanding  of  three-dimensional  separation 
and  accompanying  vortex  formation  may  help  in  adding  to  the  aerodynamic 
efficiency  of  these  aircraft.  Design  efforts  which  try  to  take  advantage  of 
vortex  lift  or  utili2e  vortex  flaps  will  surely  benefit  from  basic  studies  of 
separated  flows.  Controllability  problems  associated  with  vortex  interactions 
and  asymmetric  forebody  separation  will  benefit  from  an  understanding  of  flow 
separation.  If  separated  flows  can  be  studied  accurately  and  efficiently 
using  computational  techniques,  the  development  costs  of  these  aircraft  will 
be  significantly  reduced. 

This  report  presents  a  numerical  study  of  separation  and  accompanying 
vortices  in  the  flow  about  delta  wings  and  simple  forebodies  at  supersonic 
speeds.  The  governing  equations  considered  in  this  investigation  are  inviscid 
(i.e.,  Euler's  equation).  It  is  well  known  that  separation  zones  (i.e., 
spirals  or  closed  separation  bubbles  in  2-D  flow)  are  associated  with  the 
existence  of  vorticity  in  the  flow.^»^’^»^  There  are  two  sources  of  vorticity 
in  the  high  Reynolds-number  flows  about  supersonic  aircraft:  first,  the 
vorticity  in  the  boundary  layer,  and  second,  the  vorticity  produced  by  shock 
waves.  When  boundary  layer  separation  occurs,  the  vorticity  of  the  boundary 
layer  is  shed  into  the  flow  field  and  may,  to  a  good  approximation,  be 


considered  confined  to  an  infinitesimal  sheet  which  rolls  into  a  vortex.  This 
approximation  has  been  used  in  conjunction  with  linearized  flow  field  models 
for  many  years  (see,  for  example.  Ref  5  and  6)  v'ith  good  results.  The  present 
author^  showed  preliminary  Euler  calculations  which  indicated  that  the  more 
exact  flow  field  model  made  the  computation  of  vortex  flows  simpler  and  at  the 
same  time  more  accurate.  At  flight  Reynolds  numbers  (>  10®),  viscous  effects 
are  confined  to  thin  flow  regimes  (boundary  layers  and  vortex  sheets)  with  the 
majority  of  the  flow  inviscid  in  nature.  This  fact  makes  the  solution  of  the 
Navier-Stokes  equations  with  a  unified  procedure  (no  special  treatment  of 
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sheets  or  boundary  layers)  very  difficult.  The  Euler  equations  combined  with 
special  treatments  for  sheets  and  boundary  layer  may  turn  out  to  be  a  more 
reliable  and  accurate  approach  to  the  prediction  of  these  flows.  The  power  of 
the  Euler  equation  model  lies  In  its  ability  to  transport  vorticity  while 
avoiding  the  evaluation  of  the  boundary  layer  details.  The  extent  to  which 
the  Euler  equations  can  be  used  to  predict  separated  flows  Is  a  major  goal  of 
this  work. 

In  the  past  several  years,  a  number  of  Investigators  have  studied 
numerical  solutions  to  the  Euler  equations  exhibiting  separation.  In  Ref  8, 
Salas  showed  how  shock  vorticity  can  cause  separation  at  the  base  of  a 
cylinder.  Rizzi  et  al^  and  Murman^®  have  shown  solutions  to  Euler's  equations 
with  leading  edge  separation.  In  Ref  9  and  10,  the  source  of  vorticity  is 
unclear;  there  is  no  evidence  of  strong  shock  waves.  Both  authors  indicate 
the  possibility  that  separation  is  caused  by  numerical  viscosity.  In  Ref  10, 
the  wing  is  a  flat  plate  delta;  therefore,  there  must  be  a  singularity  at  the 
leading  edge.  The  numerical  scheme  used  could  implicitly  introduce  a  Kutta 
condition  at  the  leading  edge  and  shed  vorticity  into  the  flow  field.  R. 
Newsome^^  showed  Euler  solutions  over  an  elliptic  delta  wing  with  finite 
thickness.  His  crude  grid  results  exhibited  leading  edge  separation  (before 
any  crossflow  shock);  his  very  fine  grid  results  exhibited  separation  after  a 
crossflow  shock.  The  numerical  schemes  of  Ref  9,  10,  and  11  require  the 
addition  of  artificial  viscosity  for  stability.  Their  research  has  indicated 
a  small  sensitivity  of  separation  zones  to  reductions  in  artificial 
viscosity.  It  seems  that  these  separation  zones  are  Reynolds  number 
independent  (even  numerical  Reynolds  number  independent)  and  one  must  be  very 
careful  in  controlling  artificial  viscosity.  It  is  generally  agreed  today 
that  the  Euler  separation  of  the  type  shown  in  Ref  9,  10,  and  11  is  due  to 
explicitly  added  dissipation.  It  is  Interesting  to  note  that  the  separated 
flows  of  Ref  9,  10,  and  11  compare  well  (at  least  qualitatively)  with 
experimental  data  even  though  the  source  of  vorticity  is  numerical  in  nature, 
particularly  when  the  separation  point  is  obvious,  as  in  the  case  of  the  flat- 
plate  delta  wings  of  Ref  10. 

In  the  Euler  calculations  presented  here,  separation  occurs  only  when  a 
well-defined  source  of  vorticity  is  present.  In  addition  to  vorticity  shed 
from  the  boundary  layer,  shock  waves  in  supersonic  conical  flow  can  produce 


enough  vorticity  to  induce  separation.  This  type  of  separation  is 
qualitatively  similar  to  that  produced  by  vorticity  shed  from  a  boundary 
layer.  With  boundary  layer  vorticity  excluded,  any  comparison  with 
experimental  data  will  be  poor.  It  is  difficult  to  assess  the  relative 
magnitude  of  the  vorticity  introduced  into  the  flow  from  shocks  and  the 
boundary  layer.  It  should  be  pointed  out  that  experimental  data  show  large 
regions  of  separated  flow  with  no  significant  shock  vorticity.  This  seems  to 
indicate  that  vorticity  shed  from  the  boundary  layer  is  much  larger  than  that 
produced  by  shocks.  In  addition,  boundary  layer  separation  has  a  tendency  to 
reduce  shock  strengths  and,  therefore,  shock  vorticity.  Nevertheless,  the 
calculations  presented  here  show  that  shock  vorticity  may  play  an  important 
role  in  the  separation  process.  The  author  investigated  the  effects  of  shock 
vorticity  on  circular  cones^*^^  (Ref  12  is  included  in  Appendix  A  for 
completeness)  and  on  elliptic  delta  wings, and  found  that  at  high  angle  of 
attack  a  crossflow  shock  (Fig.  1)  can  produce  enough  vorticity  to  cause 
separation.  The  interaction  of  the  shock  and  vortex  system  will  be  discussed 
in  detail  later  in  this  report.  In  addition,  the  relationship  between  shed 
and  shock  vorticity  will  be  outlined. 

Most  geometries  considered  in  this  work  are  either  delta  wings  with 
similar  cross  .sections  or  circular  cones,  so  that  each  section  normal  to  the  z 
axis  (Fig.  1)  is  self  similar  and  the  geometry  is  conical.  The  flow  is 
assumed  supersonic  everywhere  and  separation  lines  are  assumed  to  be  on 
conical  rays.  With  these  assumptions,  it  follows  that  the  flow  is  conical 
(i.e.,  all  flow  variables  are  independent  of  the  spherical  radius  r  (Fig. 

1)).  However,  the  basic  features  of  the  flow,  including  separatior>  and 
vortices,  are  three-dimensional  in  nature.  In  addition,  the  models  used  here 
have  been  extended  to  fully  three-dimensional  flow  and  these  results  will  be 
presented  later  in  this  report. 

Our  previous  work  in  highly  vortical  flows^’^^’^^  will  not  be  presented 
in  this  report,  as  this  was  covered  before  in  Part  I.  The  overall 
computational  procedure  will  be  outlined  and  the  model  for  shedding  vorticity 
from  specified  separation  points  will  be  discussed.  Flow  fields  with  shock- 
vorticity-induced  separation  will  be  studied,  in  addition  to  flows  with 
vorticity  shed  from  both  primary  and  secondary  boundary  layer  separation 
points.  A  detailed  comparison  of  these  two  different  sources  of  vorticity 


will  be  presented.  The  investigation  of  the  computed  flow  fields  will  be 
aided  by  comparisons  with  experimental  data  and  computations  of  other 
researchers.  Finally,  the  findings  of  this  work  will  be  summarized. 


2.  COMPUTATIONAL  PROCEDURE 


2.1  OVERALL  COMPUTATIONAL  PROCEDURE,  SHOCK  VORTICITY  ALONE 

Although  employing  a  number  of  new  features,  the  overall  numerical 
procedure  used  in  this  study  is  essentially  standard. The  fully  three- 
dimensional  Euler  equations  are  solved  with  an  explicit  marching  technique. 

The  marching  direction,  z  (Fig.  1),  is  an  iterative  coordinate  for  the  conical 
flow.  The  scheme  is  restricted  by  the  fact  that  the  axial  component  of 
velocity,  W  (Fig.  1),  must  be  supersonic  everywhere.  The  marching  scheme  is 
continued  until  the  flow  field  is  invariant  with  respect  to  the  computational 
marching  direction,  except  for  a  scale  factor  for  conical  flow,  and  to  the  end 
of  the  body  considered  for  fully  three-dimensional  flows.  The  finite 


difference  scheme  used  is  Moretti's  characteristic-based  \-scheme. 


The  bow 


shock  and  the  primary  crossflow  shock  are  fit  and  are  forced  to  satisfy  the 
exact  Rankine-Hugoniot  jump  conditions.  The  bow  shock  is  fit  as  the  outermost 
boundary  of  the  flow  field.  On  the  low-pressure  side  of  the  bow  shock, 
freestream  conditions  exist.  The  crossflow  shock  is  fit  as  an  internal 
boundary  of  the  flow  field  with  its  low-pressure  side  being  computed  as  the 
computation  proceeds.^® 

The  crossflow  shock  computation  is  a  critical  part  of  the  overall 
procedure,  particularly  since  it  plays  such  a  critical  role  in  separation.  As 
mentioned  in  Section  1,  it  is  the  vorticity  produced  by  the  crossflow  shock 
that  is  significant.  Any  scheme  that  captures  the  crossflow  shock  and 
introduces  additional  artificial  viscosity  to  stabilize  it  runs  the  risk  of 
distorting  the  separation  and  vortex.  The  bow  shock  also  introduces  vorticity 
into  the  flow  field,  but  not  enough  to  produce  separation.  The  crossflow 
shock  is  caused  by  the  fact  that  at  a  large  incidence  the  crossflow  component 
of  the  velocity,  /  u^+v^  (Fig.  1),  becomes  supercritical.  The  flow  expands 
in  going  around  the  body,  and  the  crossflow  can  become  supercritical  if  the 
incidence  is  high  enough.  This  component  of  velocity  must  be  small  in  the  lee 
symmetry  plane  due  to  the  boundary  conditions.  A  supercritical  crossflow 
generally  passes  through  a  relatively  normal  shock  before  stagnating  in  the 
lee  plane  on  the  surface.  The  crossflow  shock  exists  to  turn  the  three- 
dimensional  flow  parallel  to  the  lee  symmetry  plane.  The  variation  in 
strength  of  the  crossflow  shock  can  be  quite  large.  The  shock  is  strongest  at 
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the  body  where  the  crossflow  is  highly  supersonic  and  approaches  zero  strength 
in  the  field  where  the  crossflow  becomes  sonic.  It  is  this  variation  in 
strength  that  produces  a  crossflow  entropy  gradient  and  thus  radial 
vorticity.  This  vorticity  in  turn  causes  shock- induced  separation. 

The  low-pressure  side  of  the  crossflow  shock  is  computed  with  one-sided 
differences  away  from  the  shock.  This  is  consistent  with  the  fact  that  the 
crossflow  is  supersonic.  The  Rankine-Hugoniot  conditions,  together  with  the 
compatibility  condition  along  a  bicharacteristic  reaching  the  shock  on  its 
high-pressure  side,  supply  enough  information  to  compute  the  deviation  of  each 
shock  point  from  a  conical  ray  in  addition  to  all  the  primitive  variables  on 
its  high-pressure  side.  The  bicharacteristic  used  is  the  one  in  the  plane 
containing  the  local  normal  to  the  shock  and  the  marching  direction.  All 
shock  points  are  computed  with  the  post-correction  scheme  proposed  by  Rudman^^ 
and  independently  by  deNeef.^®  The  shocks  converge  (i.e.,  each  shock  point 
becomes  aligned  with  a  conical  ray)  with  the  rest  of  the  flow  field.  In  all 
the  computations  presented  here,  the  last  crossflow  shock  point  fit  had  a 
normal  Mach  number  of  approximately  1.05  (pressure  ratio  of  1.12).  The  finite 
difference  scheme  was  able  to  capture  weaker  shock  points.  In  addition,  the 
finite  difference  scheme  was  used  to  capture  all  reverse  crossflow  shocks  and 
oblique  crossflow  shocks  in  the  flow  field  (to  be  discussed  later).  These 
shocks  are  usually  weak  and  the  scheme  can  capture  them  accurately. 

The  computational  grid  used  in  this  work  is  developed  in  stages.  In  the 
first  stage,  the  flat  wing  section  in  a  constant  z  plane  is  mapped  to  a 
circle.  This  step  is  omitted  in  the  case  of  a  circular  cone.  The  mapping  is 
a  simple  Joukowski  transformation  and  will  not  be  detailed  here.  Next,  a 
polar  coordinate  system  is  used  in  the  mapped  space.  An  exponential 
stretching  is  used  in  the  radial  direction  to  cluster  grid  points  near  the 
surface  in  order  to  resolve  the  complex  flow  near  separation.  The  Joukowski 
mapping  clusters  grid  points  in  the  circumferential  direction  automatically. 
When  thin  ellipses  (axis  ratios  >  10)  were  first  considered, solutions  with 
large  truncation  error  at  the  wing  leading  edge  were  found.  One  result  of 
this  truncation  error  was  the  production  of  vorticity  at  the  leading  edge. 
There  was  little  error  in  entropy  at  the  leading  edge  (i.e.,  entropy  was 
constant  on  streamlines),  but  the  vorticity  was  large,  leading  to  a  violation 
of  Crocco's  theorem.  The  error  was  manifested  by  a  layer  in  crossflow 
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velocity  at  the  surface.  The  crossflow  on  the  body  was  much  lower  than  that 
one  point  away  from  the  body.  While  this  vorticity  was  never  large  enough  to 
cause  leading  edge  separation,  as  in  the  work  of  Ref  9,  10,  and  11,  it  did 
distort  the  separation  behind  the  crossflow  shock.  This  error  was  eliminated 
by  increasing  the  resolution  at  the  wing  leading  edge.  An  exponential 
stretching  was  used  in  the  circumferential  direction  in  the  delta  wing  results 
shown  here. 

Shock-induced  inviscid  separation  produces  a  contact  surface  emanating 
from  the  separation  point  (Fig.  1).  It  is  this  contact  that  ultimately 
spirals  to  form  a  vortex.  The  contact  surface  has  a  jump  in  entropy  and  an 
accompanying  jump  in  velocity.  In  conical  flow,  entropy  is  constant  on 
crossflow  streamlines.  Thus,  one  can  see  (Fig.  1)  that  the  entropy  on  the 
windward  side  of  the  contact  comes  from  the  windward  stagnation  point  of  a 
cone  and  the  wing  leading  edge  stagnation  point  in  the  case  of  an  ellipse. 

The  streamline  that  wets  the  body  passes  through  the  base  of  the  crossflow 
shock  to  form  the  high-entropy  side  of  the  contact.  There  are  a  number  of 
possibilities  for  the  entropy  on  the  lee  side  of  the  contact.  The  crossflow 
streamlines  that  wet  the  body  are  tracked  in  each  step  of  the  iteration,  and 
the  proper  entropy  is  imposed  on  the  surface,  including  the  jump  in  entropy  at 
the  separation  point.  All  entropy  discontinuities  off  the  body  are  captured 
with  the  finite  difference  scheme. 

The  finite  difference  scheme  used  in  this  work  is  an  explicit  marching 
scheme  which  is  notoriously  inefficient  for  converging  to  a  conical  or  steady 
solution.  One  advantage  of  the  explicit  marching  scheme  is  that  it  is  totally 
" vector i z-ab  1  e. "  The  time-consuming  parts  of  the  code  (i.e.,  interior  point 
computation)  utilized  the  vector  architecture  of  the  Cray  1  computer.  This 
made  the  computation  about  20  times  faster  than  on  an  IBM  3033.  The 
computations  shown  here  typically  took  20  CPU  minutes  on  the  Cray  XMP,  and  the 
conical  grid  used  had  89  x  89  points  in  the  cross  sectional  plane.  In  each 
case  shown,  the  maximum  residual  (i.e.,  the  derivative  of  pressure  in  the 
marching  direction)  was  reduced  at  least  six  orders  of  magnitude. 

2.2  SHED  VORTICITY  -  COMPUTATIONAL  MODELS 

In  this  section,  the  concept  of  shedding  vorticity  from  a  specified 
separation  point  (the  intersection  of  a  conical  separation  line  and  the  cross 
plane)  will  be  discussed.  The  conceptual  process  of  shedding  vorticity  from  a 


conical  surface  will  force  separation  and  an  accompanying  spiral.  The  model 
used  here  to  force  separation  at  a  specified  location  is  founded  in  the  work 
of  J.H.B.  Smith.®  Smith  assumed  irrotational  flow  outside  a  vortex  sheet  in 
order  to  analyze  the  local  flow  at  separation.  The  present  work  uses  only  the 
basic  concept,  which  does  not  depend  on  the  irrotational  assumption.  In 
addition,  any  model  of  inviscid  separation  should  reduce  to  that  of  Smith  as 
the  rotational ity  outside  the  vortex  sheet  gets  small.  The  basic  concept 
acquired  from  Ref  5  is  that,  at  a  specified  separation  point,  there  is  a 
vortex  sheet  or  contact  surface  which  has  a  jump  in  velocity  direction.  The 
condition  at  the  intersection  of  the  sheet  and  a  smooth  surface  (i.e.,  the 
separation  point)  is  that  the  crossflow  velocity  stagnates  on  the  lee  side  of 
the  sheet  (Fig.  2a).  The  concept  of  crossflow  applies  only  to  conical  flow. 
The  extension  of  this  model  to  fully  3-0  flow  will  be  discussed  in  another 
part  of  this  section.  The  crossflow  velocity  on  the  wind  side  of  the  sheet  is 
determined  by  the  global  solution.  In  general,  the  crossflow  velocity  on  the 
wind  side  of  the  sheet  is  finite.  The  vortex  sheet  is  a  stream  surface  and 
the  pressure  is  continuous  across  it.  In  the  case  of  isentropic  irrotational 
flow  considered  in  Ref  5,  these  conditions  imply  that  the  modulus  of  velocity 
is  continuous  across  the  sheet.  If  the  crossflow  stagnates  on  both  sides  of 
the  sheet  (u  =  v  *  0)  at  separation,  the  isentropic  condition  implies  that  the 
radial  component  of  velocity  (w.  Fig.  1)  is  continuous  and  the  sheet  does  not 
exist,  i.e.,  no  vorticity  is  shed  into  the  flow  field.  These  arguments 
conclude  that,  in  isentropic  flow,  the  separating  sheet  must  leave  the  cone 
surface  tangentially  if  the  wing  side  crossflow  is  subsonic.  If  the  surface 
crossflow  is  subsonic  just  upstream  of  separation  and  the  sheet  leaves  the 
surface  at  any  angle  relative  to  the  surface,  then  the  flow  on  the  windside  of 
the  sheet  must  stagnate  and  no  vorticity  will  be  shed  into  the  flow.  On  the 
other  hand,  if  the  crossflow  is  supersonic,  the  sheet  can  leave  the  surface  at 
an  angle,  an  oblique  crossflow  shock  will  occur,  and  the  flow  on  the  wind  side 
of  sheet  will  not  stagnate  so  that  vorticity  is  shed.  This  phenomenon  will  be 
discussed  further  in  this  section.  None  of  these  arguments  holds  for  the 
highly  nonisentropic  flow  discussed  in  the  previous  section,  where  the  flow 
stagnates  on  both  sides  of  the  separating  sheet.  In  these  cases,  vorticity  is 
not  shed  from  the  surface;  only  shock  vorticity  is  present. 

Separation  was  forced  in  the  calculations  presented  here  simply  by  using 
a  double  point  at  separation  in  order  to  allow  for  a  jump  in  crossflow 


velocity.  The  crossflow  velocity  (u)  on  the  wind  side  of  the  sheet  is 
determined  by  the  governing  crossflow  momentum  equation,  as  at  any  other 
surface  point.  The  crossflow  component  (v)  must  vanish  at  every  surface  point 
from  the  boundary  condition.  The  crossflow  velocity  on  the  lee  side  of  the 
sheet  is  set  to  zero  (u  =  0)  in  order  to  force  separation  at  a  specified 
point. 

The  separation  point  in  this  forced  separation  model  exhibits  a  Jump  in 
crossflow  velocity.  The  crossflow  on  the  wind  side  (Fig.  23a)  of  the 
separation  is  evaluated  in  each  step  of  the  iteration  from  the  momentum 
equation  tangent  to  the  body  in  the  cross  plane.  The  remaining  primitive 
variables  are  computed  at  this  point  with  the  same  governing  equations  as  at 
any  other  body  point,  with  the  boundary  condition  (v  »  0)  satisfied.  All 
circumferential  velocity  derivatives  are  taken  one-sided  in  the  negative  9 
direction  in  order  to  avoid  differencing  across  the  sheet.  The  point  just 
below  the  sheet  (lee  side)  is  forced  to  be  a  crossflow  stagnation  point 
(u  =  V  =  0).  The  pressure  is  continuous  across  the  sheet,  so  that  its  value 
on  the  wind  side  could  theoretically  be  used  on  the  lee  side  at  separation. 
Difficulties  were  encountered  in  computing  the  pressure  in  a  small  region  just 
after  separation.  The  computation  of  the  pressure  in  this  region  will  be 
discussed  in  the  next  paragraph.  The  entropy  at  the  lee  side  point  is 
computed  in  the  standard  way,  i.e.,  conserved  along  the  streamline  that  wets 
the  lee  side  of  the  body.  If  the  pressure  is  known  in  addition  to  the  entropy 
and  total  temperature  (the  flow  is  assumed  adiabatic),  the  last  component  of 
velocity  (w)  (Fig.  1)  can  be  evaluated  from  the  energy  equation.  Again, 
circumferential  derivatives  of  the  velocities  across  the  sheet  at  the  body 
grid  points  near  separation  are  avoided.  Differences  in  entropy  across  the 
sheet  are  avoided  naturally  by  the  x-scheme,  as  in  the  shock- induced 
separation  discussed  in  the  previous  section. 

The  difficulty  encountered  in  the  evaluation  of  the  pressure  on  the  lee 
side  of  the  separation  point  can  be  traced  to  the  numerical  computation  of  the 
derivatives  of  velocities  in  a  direction  normal  to  the  surface.  The 
derivatives  appear  in  the  continuity  equation  from  which  the  pressure  is 
evaluated.  These  velocities  are  discontinuous  across  the  sheet  and,  since  the 
sheet  lies  so  close  to  the  surface  at  separation,  there  is  no  way  to  resolve 
the  flow  in  the  region  between  the  surface  and  the  sheet  near  separation.  The 
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behavior  of  the  pressure  1n  this  region  was  evaluated  analytically  by  Smith- 


and  numerically  by  Fiddes®  under  the  assumption  of  slender-body  theory.  The 
results  of  Fiddes  show  a  pressure  plateau  on  the  surface  just  after 
separation.  The  full  Euler  equations  predict  that  the  pressure  derivative 
along  the  surface  (pg)  becomes  zero  at  the  lee  side  of  the  separation  point, 
where  the  crossflow  stagnates.  In  the  present  calculation  it  has  been 
difficult  to  obtain  this  type  of  pressure  behavior.  A  number  of  procedures 
have  been  attempted  with  varying  degrees  of  success,  but  none  has  been  totally 
satisfactory.  In  Ref  7  the  author  showed  preliminary  results  which  exhibited 
a  pressure  plateau  after  separation,  but  the  scheme  used  there  proved  unstable 
when  the  vortex  sheet  was  moved  during  the  Iteration  process.  A  reformulation 
of  dependent  variables  has  alleviated  the  problem.  Currently,  as  long  as  the 
crossflow  Is  subsonic  (Fig.  2a)  just  before  separation,  no  special  treatment 
is  used  in  evaluating  the  pressure  except  for  circumferential  and  radial  grid 
clusterings.  Because  of  the  difficulties  associated  with  the  pressure 
evaluation  in  the  region  on  the  surface  just  after  separation,  the  results  of 
the  next  section  should  be  considered  somewhat  preliminary. 


In  the  case  of  supersonic  crossflow  (Fig.  2b),  differencing  of  pressure 
across  the  separation  point  proved  to  be  unstable.  The  results  of  a  number  of 
numerical  experiments  indicate  that  the  model  for  the  inviscid  flow  in  this 
case  should  be  such  that  the  sheet  leaves  the  surface  at  an  angle  relative  to 
it  and  an  associated  shock  occurs  at  separation.  The  flow  structure  is 
sketched  in  Fig.  2b.  This  model  does  not  violate  the  concepts  of  vorticity 
shedding  proposed  by  Smith,®  since  the  flow  on  the  wind  side  of  the  separation 
point  need  not  stagnate,  and  thus  a  jump  In  velocity  exists  across  the 
sheet.  In  order  to  capture  the  shock  at  separation  in  one  mesh  interval,  no 
pressure  derivatives  were  taken  across  the  separation  point  on  the  cone 
surface.  In  addition  a  small  pressure  plateau  is  imposed  just  after 
separation,  the  level  being  taken  from  just  downstream  of  separation. 

Secondary  separation  is  forced  by  shedding  vorticity  from  the  body  in  the 
reverse  crossflow  region  produced  by  primary  separation.  The  model  is  the 
same  as  was  just  outlined  for  primary  separation.  Both  separation  point 
locations  are  obtained  from  experimental  data.  The  only  free  parameters  in 
the  problem  as  posed  here  are  the  locations  of  the  separation  points.  In  the 
work  of  Ref  6,  primary  separation  points  were  found  Iteratively  by  matching  an 
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inviscid  solution  (slender-body  theory)  to  a  boundary  layer  solution.  The 
procedure  just  outlined  Is  directly  applicable  to  fully  three-dimensional 
flows  If  the  crossflow  Is  replaced  by  the  flow  In  a  plane  locally  normal  to 
the  separation  line.  Consider  Fig.  3,  where  the  separation  point  location  in 
each  plane  Z  =  constant  Is  specified  (I.e.,  experimentally).  Given  all  these 
separation  locations,  a  separation  line  can  be  constructed  (see  sketch  in  Fig. 
3).  In  each  marching  step  of  the  calculation,  the  plane  normal  to  separation 
line  can  be  determined  and  the  model  just  outlined  can  be  applied  locally  in 
each  of  these  planes.  The  zero-crossflow  condition  becomes  equivalent  to 
forcing  the  surface  velocity  vector  parallel  to  the  separation  line,  as  it 
should. 

2.3  SEPARATION  CONDITION  AT  SHARP  WING  LEADING  EDGE 

The  flow  at  the  leading  edge  of  the  zero-thickness  delta  wing  is 
singular.  Powell  and  Murman  have  been  successful  In  computing  these  flows 
with  no  consideration  of  this  singularity.^^ 

Salas^^  investigated  the  conical  flow  with  a  singular  expansion  corner. 
The  flow  at  the  wing  leading  edge  Is  similar.  Figure  4  depicts  the  flow 
situation  1n  the  crossflow  plane.  The  singularity  is  two-dimensional  in  a 
plane  normal  to  the  wing  leading  edge  (the  crossflow  plane)  since  all 
derivatives  normal  to  this  plane  are  zero  for  a  conical  flow.  The  velocity 
component  normal  to  this  plane  is  unchanged  through  the  singularity. 

The  crossflow  stagnates  just  inboard  of  the  leading  edge  on  the  wind  side 
of  the  sheet.  It  rapidly  accelerates  to  sonic  at  the  leading  edge.  The  sonic 
line  is  followed  by  a  centered  Prandtl-Meyer  expansion.  The  flow  cannot 
expand  through  a  deflection  5  >  =  130.5  .  At  a  deflection  of  130.5°, 

the  Mach  number  would  be  infinite  and  the  pressure  zero.  Hence,  the  flow 
separates  to  form  a  vortex  sheet  that  spirals  downstream.  The  deflection 
angle  6  is  such  that  the  pressure  on  the  wind  side  of  the  sheet  matches  that 
at  the  crossflow  stagnation  point  on  the  lee  side  of  the  sheet. 

During  the  computation,  a  is  evaluated  from  the  local  flow  directions, 
expanding  the  sonic  crossflow  through  the  Prandtl-Meyer  to  get  the  pressure  on 
either  side  of  the  sheet  at  the  leading  edge.  The  singular  behavior  is 
imposed  at  each  step  of  the  iteration,  and  the  deflection  6  is  converged  with 
the  rest  of  the  flow  field. 
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Fig.  3  Sketch  of  Three-dimensional  Separated  Flow 


3.  COMPARISON  OF  SHOCK  VORTICITY  AND  SHED  VORTICITY 


In  an  effort  to  gain  a  better  understanding  of  the  relationship  between 
!  shock  vorticity  and  vorticity  shed  from  the  surface  of  a  smooth  body,  the  flow 

about  the  5“  cone  tested  by  Rainbird^^  was  computed  with  a  number  of  different 
I  separation  point  locations  specified.  The  case  considered  was  =  4.25  and 

a  =  12.35°.  This  flow  is  supercritical  and,  with  no  vorticity  shed  from  the 
cone  surface,  the  crossflow  shock  produces  enough  vorticity  to  cause 
separation.  Figures  5  and  6  show  the  crossflow  streamlines  and  isobars, 
respectively,  for  this  flow  with  no  vorticity  shed  from  the  body.  Separation 
for  this  case  is  computed  to  be  at  e  =  151.3°.  Figure  5  indicates  that  the 
separating  streamlines  leave  the  surface  at  a  large  angle  (57°)  relative  to 
it.  As  shown  by  the  author  in  Ref  12,  when  only  shock  vorticity  is  present, 
there  is  no  jump  in  crossflow  velocity  at  the  separation  point,  which  is 

consistent  with  the  fact  that  no  vorticity  is  being  shed  from  the  surface.  It 

should  be  pointed  out  that  in  the  computational  results  in  this  section  all 
crossflow  shocks  are  captured.  Figure  6  indicates  that  the  shock  is  captured 
very  sharply  (see  the  closely  spaced  isobars).  Additionally,  these  captured 
shock  results  compare  very  well  with  the  shock  fit  results  for  this  case  (Ref 
12). 

Figures  7  and  8  show  results  for  the  other  extremes  of  separation  point 
location  studied  (9  =  115°).  Figure  7  shows  the  crossflow  streamlines.  The 
secondary  separation  shown  is  due  to  a  strong  reverse  crossflow  shock  (see 

Fig.  8),  and  the  third  vortex  off  the  surface  is  similar  to  the  one  discussed 

previously.  In  Fig.  8,  the  isobars  are  shown,  and  they  indicate  an  oblique 
shock  at  the  specified  primary  separation  point.  The  jump  in  velocity  at  the 
separation  point  is  significant  with  separation  specified  at  9  =  115°, 
indicating  significant  vorticity  being  shed  from  the  surface.  A  comparison  of 
Fig.  5  and  7  shows  that  the  extent  of  the  vortical  regions  is  comparable, 
while  the  two  sources  of  vorticity  are  very  different. 

The  relationship  between  shock  vorticity  and  shed  vorticity  is  made 
clearer  by  considering  Fig.  9.  The  figure  shows  the  jump  in  crossflow 
velocity  vs  separation  point  location.  The  jump  in  crossflow  velocity  is 
directly  related  to  the  vorticity  shed  into  the  flow  field  from  the  separation 
point.  The  shock  configuration  transition  from  an  oblique  crossflow  shock  to 


a  detached  normal  crossflow  shock  occurs  at  about  a  128”  (indicated  by  the 
shaded  area  in  Fig.  9).  It  should  be  pointed  out  that  the  jump  in  velocity  at 
the  separation  point  in  the  oblique  shock  cases  was  computed  by  subtracting 
the  oblique  shock  velocity  jump  from  the  numerical  results.  Thus,  the  jumps 
in  velocity  in  Fig.  9  represent  the  jumps  across  the  vortex  sheet  at 
separation.  The  figure  shows  that  this  velocity  jump  goes  to  zero  smoothly  as 
the  separation  point  location  due  to  shock  vorticity  alone  is  approached  (e^  = 
151.3°).  This  indicates  that  separation  due  to  shock  vorticity  alone  and  that 
due  to  shed  vorticity  are  related.  In  fact,  it  would  seem  that  separation  due 
to  shock  vorticity  alone  is  a  particular  solution  of  the  set  of  solutions  in 
which  vorticity  is  shed  from  the  surface.  In  this  particular  solution,  the 
value  of  the  vorticity  shed  is  zero.  Additional  details  of  this  investigation 
can  be  found  in  Ref  21,  which  is  included  here  as  Appendix  B. 
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4.  COMPLEX  SHOCK  PATTERNS  AND  VORTICES 


In  an  investigation  of  high  Mach-nuinber,  high  angle-of-attack  flows,  a 
number  of  interesting  flow  phenomena  were  discovered.  These  were  due  to  the 
very  high  vorticity  imparted  into  the  flow  from  the  shock  system  and  the 
interaction  of  the  resulting  vortices  with  that  same  shock  system. 

In  this  computational  investigation,  a  single  cone  is  considered.  The 
half  angle  is  10°,  and  two  Mach  numbers  (5  and  10)  were  studied.  At  a 
freestream  Mach  number  of  5,  the  shock  system  becomes  complex  at  an  angle  of 
attack  of  25°.  Figure  10  shows  the  crossflow  streamlines  and  Fig.  11  the 
isobars.  The  crossflow  shock  exhibits  a  slope  discontinuity  as  it  passes  over 
the  top  of  the  vortex,  which  it  created.  This  discontinuity  in  shock  slope  is 
accompanied  by  an  additional  shock,  which  extends  toward  the  lee  plane.  This 
shock  is  captured  but  is  clearly  seen  in  the  isobar  plot  of  Fig.  11.  The 
crossflow  shock  is  no  longer  simple,  but  has  a  triple  point  separating  two 
portions  of  shock. 

Around  the  primary  vortex  a  ''pinwheel'*^^  shock  pattern  is  formed.  The 
crossflow  shock  becomes  oblique  to  the  flow  just  off  the  surface  so  that  the 
crossflow  behind  the  shock  is  supersonic  just  off  the  surface  (see  the  sonic 
line  plot  of  Fig.  12).  The  bundle  of  supersonic  crossflow  extends  to  the 
triple  point.  This  flow  must  be  shocked  in  order  to  meet  the  symmetry 
condition  in  the  lee  plane,  and  so  a  crossflow  shock  sits  on  top  of  the  vortex 
(Fig.  11).  As  lee  plane  flow  is  accelerated  by  the  vortex,  it  becomes 
supersonic,  so  that  another  sonic  bubble  is  formed  attached  to  the  lee  plane 
(Fig.  12).  This  supersonic  flow  must  be  shocked  to  meet  the  surface  boundary 
condition.  Note  in  Fig.  11  that  the  shock  is  normal  to  the  lee  plane  just 
above  the  body.  Finally,  the  surface  flow  is  accelerated  under  the  primary 
vortex  and  becomes  supersonic  (see  the  third  sonic  bubble  in  Fig.  12).  This 
flow  is  shocked  in  order  to  become  tangent  to  the  inside  of  the  separation 
line.  This  process  of  expanding  and  shocking  as  the  flow  moves  around  the 
vortex  sets  up  a  system  of  three  shocks  which  the  author  first  referred  to  as 
the  pinwheel  type  in  Ref  22.  This  shock  system  is  clearly  shown  in  the 
isobars  of  Fig.  11. 
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Complex  shock  patterns  occur  at  a  =  20“  for  =  10  instead  of  a  =  25“  at 
^co  ■  5.  Figures  13  and  14  show  the  crossflow  streamlines  and  isobars, 
respectively,  for  the  =  10  and  a  =  20“  case.  Figure  14  shows  the  pinwheel- 
type  shock  pattern  discussed  previously,  in  addition  to  the  triple  point  in 
the  crossflow  shock.  It  should  be  noted  that  there  is  no  secondary  vortex  in 
this  case;  instead,  there  is  a  simple  node  in  the  lee  plane  above  the  primary 
vortex. 

The  last  case  considered  is  =  10  and  a  =  25°.  Figure  15  shows  the 
streamlines  and  Fig.  16  the  isobars.  Again,  the  pinwheel  shock  pattern  is 
quite  evident  in  addition  to  the  triple  point.  The  primary  crossflow  shock 
has  a  maximum  normal  Mach  number  of  3.8  at  the  surface  and  is  quite  strong. 

The  captured  reverse  crossflow  shock  has  a  normal  Mach  number  of  1.8.  This 
has  a  secondary  vortex  above  the  main  vortex,  in  addition  to  a  simple  node  in 
the  lee  plane  halfway  between  the  body  and  the  bow  shock  (not  shown  in  the 
figure).  The  secondary  spiral  node  above  the  main  vortex  is  seen  more  clearly 
in  the  expanded  view  of  its  region  shown  in  Fig.  17.  Figure  17  is  an  expanded 
view  of  the  dashed  box  in  Fig.  15.  It  can  be  seen  that  the  node  does  not  form 
a  large  vortex  as  does  the  primary  vortex,  but  seems  to  spiral  at  the  node. 
Additional  details  of  this  work  can  be  found  in  Ref  23,  which  is  included  here 
as  Appendix  C. 
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5.  DELTA  WING  FLOWS 


Over  the  past  several  years  a  number  of  researchers  investigating 
separated  flow  about  sharp  leading  edge  wings  have  concluded  that  there  exists 
a  close  relationship  between  flow  separation  and  computed  total  pressure  loss 
(see,  for  example.  Ref  24).  This  is  not  surprising  in  light  of  Crocco's 
equation  which  relates  vorticity  to  the  gradient  of  entropy  (i.e.,  total 
pressure).  But  this  conclusion  seems  inconsistent  with  the  classical  model 
for  inviscid  separated  flow.  The  classical  model^*^*®  is  one  of  isentropic 
irrotational  flow  with  a  infinitesimally  thin  region  of  vorticity,  a  vortex 
sheet.  This  apparent  inconsistency  was  investigated,  and  it  was  found  that  a 
zero  total  pressure  loss  Euler  solution  with  leading  edge  separation  does 
exist. 

In  all  the  cases  considered  in  this  investigation,  the  shock  system 
produces  very  little  entropy.  The  question  then  arises  as  to  how  Crocco's 
relation  vH^-TvS-?  x  2=0  is  satisfied  by  a  separated  flow  with  vorticity  and 
no  entropy  or  total  enthalpy  gradients.  The  exact  inviscid  solution  is 
irrotational  (2  =  0)  everywhere  outside  the  vortex  sheet  and  away  from  the 
vortex  core.  These  flows  have  been  computed  for  many  years  using  irrotational 
models  (see  Ref  25,  for  example)  with  the  sheet  and  the  vortex  core  being 
modeled  by  singularities.  Modern  Euler  solvers  capcure  the  sheet  and  the 
vortex  core  as  part  of  the  solution.  Powell  and  Murman^®  showed  two  different 
ways  that  Crocco's  relation  could  be  satisfied  by  an  Euler  solver:  a  strictly 
conservative  formulation  in  which  entropy  gradients  generated  at  the  wing 
leading  edge  are  balanced  by  7  x  2  (Hq  is  constant)  and  a  modified 
conservative  scheme  which  substituted  S  =  constant  for  the  axial  momentum 
equation.  In  the  latter  formulation,  vS  =  7  x  2  =  0  and  again  vH^  =  0.  The 
present  formulation  is  closer  to  the  second  of  these.  Since  entropy  is 
tracked  very  accurately  along  streamlines,  the  present  formulation  differs 
from  the  second  of  Powell  and  Murman's  in  that  all  momentum  equations  are 
satisfied. 

If  we  consider  a  75°  sweep  delta  at  a  low  supersonic  Mach  number  (M^  = 
1.5)  and  small  angle  of  attack  (o  =  4°),  then  the  shock  system  (only  a  bow 
shock)  produces  less  than  1%  total  pressure  loss.  The  present  calculation 
predicts  that  result  very  accurately.  Figure  2a  shows  contour  lines  of  the 
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modulus  of  vorticity,  and  Fig.  2b  shows  the  crossflow  streamlines.  The  vortex 
is  formed  because  of  the  presence  of  the  vorticity  which  comes  from  the 
leading  edge.  The  flow  separates  because  an  attached  leading  edge  requires 
expansion  beyond  the  vacuum  limit,  as  discussed  in  Section  2  of  this  report. 
The  sheet  is  spread  over  a  few  grid  cells  near  the  leading  edge  and  gets 
thicker  as  it  wraps  around  the  vortex  core.  Unlike  the  calculations  of 
Ref  24,  no  total  pressure  loss  is  associated  with  this  separation.  In  the 
core  of  the  vortex,  Crocco's  relation  is  satisfied  to  within  the  accuracy  of 
the  calculation.  The  gradient  of  total  enthalpy  is  at  a  machine  zero,  and  the 
gradient  of  entropy  vS  is  less  than  10”^^,  so  that  only  the  ^  x  ui  remains. 
This  term  has  a  local  maximum  near  the  core  of  10"^  (|ui|  =  0(1)  at  the 
core),  which  is  within  the  accuracy  of  the  calculation.  There  are  larger 
values  of  9xu(10”^)  at  the  wing  leading  edge  due  to  the  singularity  there. 
These  findings  are  similar  to  those  of  Ref  26. 

The  case  discussed  in  the  preceding  paragraph  is  one  in  which  the  flow  is 
irrotational  outside  the  vortex  sheet  and  away  from  the  vortex  core  (see  Fig. 
18a).  The  exact  solution  to  the  differential  Euler  equations  is  approached  as 
the  sheet  approaches  an  infinitesimal  region  and  the  core  approaches  a 
point.  In  the  exact  solution,  Crocco's  relation  is  not  expected  to  be 
satisfied  across  the  vortex  sheet;  velocities  are  discontinuous  and  their 
derivatives  do  not  exist.  It  can  be  shown  that  the  flow  velocity  must  be 
parallel  to  the  vorticity  vector  at  the  spiral  node  at  the  end  of  the  vortex 
sheet  (the  vortex  core).  The  exact  solution  satisfies  Crocco's  relation  by 
forcing  i?  x  u  -  0  as  the  core  is  approached.  The  vortex  core  has  the  same 
behavior  in  the  present  calculation. 

Figure  18  shows  a  comparison  of  the  surface  pressure  computed  using  the 
present  technique  and  that  computed  using  the  fully  conservative  procedure  of 
Powell  and  Murman.^^  The  wing  is  swept  75°  and  the  free  stream  conditions  are 
=  1.95  and  a  -  10°.  The  flow  is  separated  at  the  leading  edge,  as  can  be 
seen  from  the  crossflow  streamlines  of  Fig.  20a  (the  dashed  lines  are 
crossflow  stagnating  streamlines).  The  reverse  flow  becomes  supersonic  under 
the  vortex  and  is  shocked  at  x/x^^  =>  0.75;  note  the  clustering  of  isobars  in 
Fig.  20b.  The  comparison  of  Fig.  19  indicates  that  the  two  procedures  are 
predicting  the  same  static  pressure.  The  results  of  Powell  and  Murman  exhibit 
a  65%  loss  in  total  pressure  at  the  vortex  core,  while  the  shock  system 


(b)  Crossflow  Streamlines 

R87-4811-18 

Fig.  18  Vorticity  Contourt  and  Streamlines  for  Flat  Plate  Delta 
with  X  =■  75°,  =  1.7  and  a  =  4^“ 
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produces  a  maximum  total  pressure  loss  of  only  30%.  Figure  21  compares  the 
surface  crossflow  Mach  number  computed  by  the  two  calculations.  For  the  sake 
of  clarity,  this  Mach  number  is  plotted  vs  normalized  arc  length  (i.e.,  0  to  1 
is  the  wind  side  and  1  to  2  is  the  lee).  Here  again  the  comparison  is 
generally  good.  There  are  small  discrepancies  just  beyond  the  leading  edge  in 
addition  to  those  near  the  lee  plane  of  symmetry.  In  Ref  26,  Powell  and 
Murman  found  that  the  major  discrepancy  between  results  that  conserve  entropy 
and  those  that  do  not  is  in  the  radial  component  of  Mach  number.  The 
computations  presented  here  seem  to  indicate  the  same  result. 

The  flow  at  the  leading  edge  of  the  zero-thickness  delta  wing  is 
singular.  Powell  and  Murman  have  been  successful  in  computing  these  flows 
with  no  consideration  of  this  singularity.  The  results  presented  here  in  Fig. 
18-21  were  also  computed  with  no  consideration  for  this  singularity.  The  wing 
leading  edge  was  avoided,  i.e.,  no  grid  point  was  placed  there.  A  Kutta-type 
condition  was  developed  for  sharp  leading  edge  wings,  and  was  outlined  in 
Section  2  of  this  report.  Figure  22a  shows  a  comparison  of  the  surface  pres¬ 
sure  distribution  computed  with  and  without  the  singularity  included.  The 
wing  sweep  is  TS”  and  M^  =  1.7  and  a  =  12°.  The  differences  in  the  pressures 
are  very  small  and  confined  to  the  area  of  the  tip.  Figure  22b  shows  a 
comparison  of  the  computed  surface  crossflow  Mach  number  distribution.  Here 
again  the  differences  are  small  and  confined  to  the  leading  edge  region.  It 
should  be  noted  that  the  sonic  velocity  at  the  tip  (s/X|_^  =  1)  is  indicated  by 
the  symbol  (M^^p  =  1.),  as  is  the  flow  after  the  centered  expansion 
(Mqp  =  1.3).  In  all  the  results  that  follow,  the  singular  tip  solution  is 
computed  explicitly. 

In  comparing  the  results  of  inviscid  separated  flow  calculations  with 
experimental  data,  one  will  find  significant  discrepancies  under  the  main 
vortex.  A  large  part  of  these  discrepancies  is  due  to  secondary  separation. 

In  Ref  7,  this  author  presented  a  model  for  shedding  vorticity  from  a  smooth 
surface  in  order  to  force  separation  at  a  specified  location  in  an  Euler 
calculation.  In  Ref  20,  this  work  was  extended  to  include  secondary 
separation,  i.e.,  separation  of  the  reverse  flow  under  the  main  vortex. 

Figure  23  shows  a  comparison  of  computed  surface  pressure  and 
experimental  data.^^  The  wing  is  swept  75°,  M^  =  1.7  and  a  =  12°.  Also  shown 
in  the  figure  are  the  computed  results  of  Powell  and  Murman  which  were  first 
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(b)  Surface  Crossflow  Mach  Number 

R87-4811-22 

Fig.  22  Comparison  of  Results  with  and  without  Kutta  Condition  Imposed 
\  =  75°,  1.7,0=  12° 


presented  1n  Ref  30,  Both  the  Euler  computations  predict  the  lee  side 
pressures  quite  well,  except  in  the  region  under  the  main  vortex  (x/x|_^  > 

0.6).  Figure  24a  shows  the  crossflow  streamlines,  which  indicate  that  the 
vortex  is  centered  above  x/xl^  *  0.6.  Figure  24a  shows  a  separation  of  the 
reverse  flow  which  is  due  totally  to  shock  vorticity  and  is  downstream  of  the 
secondary  separation  location  noted  experimentally.  The  shock  that  causes 
this  separation  is  the  reverse  crossflow  shock  located  at  x/Xj^^  =  0.8. 

Figure  11b  shows  the  isobars  from  which  the  shock  can  be  seen  clearly. 

The  experimental  data  of  Miller  and  Wood*^^  included  cross-sectional 
vapor-screen  photographs  from  which  secondary  vortices  can  be  seen.  The 
location  of  secondary  separation  was  difficult  to  determine  from  these 
pictures,  but  an  attempt  was  made.  Figure  25  shows  a  comparison  among 
computed  surface  pressures,  secondary  separation  forced  at  x/X|_£  =  0.7,  and 
the  experimental  data.  The  figure  shows  a  significant  improvement  in  the 
comparison  in  the  pressures  for  x/X|_£  >  0.6.  From  the  comparison  of  Fig.  25, 
one  might  get  the  impression  that  the  Euler  results  with  secondary  separation 
predict  the  flow  fields  very  accurately.  Experience  with  more  detailed 
experimental  data  indicates  that  the  test  data  are  flat  for  x/X|_£  >  0.6.  The 
overexpansions  of  the  Euler  result  under  the  main  vortex  (0.6  <  x/X[_£  <  0.7) 
and  the  secondary  vortex  (0.8  <  x/x^^  <  0.9)  are  areas  where  discrepancies 
exist.  These  discrepancies  are  due  either  to  boundary  layer  thickening  or  to 
viscous  effects  not  accounted  for  by  the  boundary  layer  equations. 

Figure  26a  shows  the  crossflow  streamlines  with  secondary  separation  at 
x/x^E  =0.7.  A  comparison  of  Fig.  26a  and  26a  indicates  that  the  forced 
secondary  vortex  is  larger  than  that  caused  by  the  reverse  crossflow  shock. 

The  secondary  separation  is  forced  to  occur  in  a  supersonic  reverse  crossflow 
region.  Therefore,  a  shock  is  produced  before  the  vortex.  Figure  26b  shows 
isobars  that  indicate  the  shock  clearly  at  x/X|_£  =  0.7.  It  must  be  remembered 
that  while  this  shock's  vorticity  may  contribute  to  this  secondary  vortex,  it 
is  the  additional  vorticity  shedding  that  causes  this  secondary  separation. 
These  comparisons  are  typical;  additional  examples  can  be  found  in  Ref  31. 


.>• 

•-PT 


r 

r 


s- 

V 


44 


(a)  Crossflow  Streamlines 
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6.  FULLY  THREE-DIMENSIONAL  FLOWS 


In  this  section  the  application  of  the  vorticity-shedding  model  used  for 
conical  flows  to  fully  three-dimensional  flows  is  demonstrated.  The 
computational  procedure  has  already  been  discussed  in  Section  2  of  this 
report.  The  present  study  is  restricted  to  simple  three-dimensional 
configurations.  The  first  is  a  cone/cylinder  configuration  tested  by  E.J. 
Landrum^^  at  NASA's  Langley  Research  Center.  The  free-stream  Mach  was  2.3  and 
the  cone  half-angle  was  9.46°.  The  conical  portion  of  the  body  ended  at  Z/L  = 
.45.  The  configuration  is  shown  in  Fig.  27.  Also  shown  in  the  figure  are  the 
surface  stream  lines  for  the  20°  angle-of -attack  case.  The  primary  separation 
line  is  clearly  indicated  by  the  discontinuity  in  the  slope  of  the  streamlines 
coming  from  the  wind  side  of  the  body.  The  location  of  this  line  was  obtained 
from  the  experimental  data.  The  streamline  picture  of  Fig.  27  is  easier  to 
understand  if,  in  addition,  the  conical  crossflow  streamlines  are 
considered.  Figure  28a  shows  the  crossflow  streamlines  on  the  conical  portion 
of  the  body  (Z  <  .45).  First,  one  sees  that  the  streamlines  coming  from  the 
wind  plane  leave  the  body  at  the  separation  point  (dashed  line,  intersecting 
the  body  at  113°  from  the  wind  plane,  is  the  separating  streamline);  this 
should  be  kept  in  mind  when  studying  Fig.  27.  Under  the  main  vortex,  there  is 
a  small  secondary  vortex  which  is  caused  by  the  vorticity  produced  by  a 
reverse  crossflow  shock.  The  crossflow  isobars  are  shown  in  Fig.  28b.  This 
shock  is  indicated  by  the  clustered  isobars  under  the  main  vortex.  There  is  a 
crossflow  reattachment  line  (dashed  line.  Fig.  28a)  on  the  body  at  124°  from 
the  lee  plane.  These  features  can  be  seen  in  the  surface  streamlines  of  Fig. 
27.  This  case  and  all  those  to  follow  have  a  supersonic  crossflow  approaching 
primary  separation.  As  mentioned  in  Section  2,  an  oblique  crossflow  shock 
occurs  at  the  separation  point,  which  can  be  seen  in  the  isobars  of  Fig. 

28b.  This  oblique  shock  is  not  captured  as  sharply  as  the  normal  reverse 
crossflow  shock. 

Figure  29  shows  a  comparison  of  the  computed  surface  pressure  and  the 
experimental  data.  Figure  29a  is  for  the  conical  flow  Z/L  <  .45;  Fig.  29b  and 
29c  are  at  Z/L  =  .675  and  .925,  respectively,  of  Landrum.-"-’  The  polar  angle  e 
is  measured  from  the  wind  plane  as  shown  in  Fig.  1.  In  the  figure,  the 
primary  separation  is  denoted  by  a  jump  in  pressure  due  to  the  oblique 


a)  Crossflow  streamlines 
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Fig.  28  Crossflow  Streamlines  and  Isobars  on  Cone/Cylinder  (2/L  <  45) 
IVI.„  »  2.3,  a  =  20" 
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crossflow  shock  mentioned  previously.  The  pressure  comparison  is  good 
upstream  of  the  primary  separation;  In  addition,  the  plateau  pressure  just 
past  separation  is  predicted  quite  well.  The  lee  plane  region  is  not 
predicted  so  well.  The  pressure  in  the  lee  plane  itself  is  predicted 
accurately,  but  the  inviscid  reverse  flow  expands  rapidly  to  a  supersonic 
reverse  Mach  number  region  which  is  terminated  by  a  reverse  crossflow  shock. 
While  reverse  crossflow  shocks  have  been  noted  experimentally,^^  they  seem  too 
strong  in  inviscid  calculations.  This  is  similar  to  the  result  found  for  the 
delta  wing  calculations  discussed  in  the  previous  section. 

The  last  case  considered  here  (more  detail  can  be  found  in  Ref  22)  is  a 
missile  configuration  with  elliptical  cross  sections.  The  configuration  was 
tested  at  NASA's  Langley  Research  Center^**  and  computed  with  a  Navier-Stokes 
solver  by  Newsome.^®  The  missile  is  made  up  of  a  series  of  3:1  elliptical 
sections.  The  computation  was  performed  at  M^  =  2.5,  a  =  20°  and  8=0  (tests 
were  conducted  at  yaw).  Primary  separation  occurred  near  the  leading  edge  of 
the  configuration;  the  location  of  the  separation  line  was  obtained  from  the 
computational  results  of  Newsome.^® 

Figure  30  shows  a  comparison  of  surface  pressures  computed  with  the  Euler 
solver,  the  Navier-Stokes  solver^^  and  that  found  experimentally.^'^  Figure 
11a  shows  a  station  near  the  apex  of  the  body  (Z/L  =*  0.3).  The  separation 
point  is  indicated  by  the  compression  (shock)  at  9  »  100°  in  the  Euler 
solution.  The  comparison  before  separation  is  good  except  for  some  upstream 
influence  seen  in  both  the  experiment  and  the  Navier-Stokes  calculation  and 
absent  in  the  Euler  solution.  In  this  flow  field  as  in  the  others  presented 
here,  the  local  crossflow  is  supersonic  at  separation.  The  pressure  plateau 
after  separation  is  predicted  well  by  both  computations.  The  Euler 
computation  predicts  supersonic  reverse  crossflow  and  a  shock  at  9  »  120°. 
Surprisingly,  it  seems  that  the  Navier-Stokes  solution  has  the  same  pressure 
distribution  in  this  region.  In  this  calculation,  the  experimental  data  do 
not  seem  to  show  this  behavior,  but  its  scatter  at  9  =  135°  is  large  and  makes 
it  difficult  to  tell  whether  a  compression  does  or  does  not  occur  in  this 
region.  Figure  30  shows  stations  farther  downstream.  The  features  are 
similar  to  those  in  Fig.  29.  The  final  two  sections  shown  in  Fig.  30  show  a 
sharp  shock  at  the  primary  separation  point,  in  addition  to  a  sharp  reverse 
crossflow  shock.  At  this  station  the  Navier-Stokes  computation  shows  no  sign 


of  supersonic  reverse  crossflow.  It  should  be  pointed  out  that  secondary 
separation  was  noted  in  the  experiment  and  Navier-Stokes  calculation  was  not 
included  in  the  Euler  computation  for  the  reasons  discussed  previously.  Here 
again,  the  Euler  computation  predicts  the  separated  flow  surface  pressures 
well,  except  in  the  region  under  the  vortex. 


7.  SUBSONIC  FLOWS 


The  separation  model  developed  here  was  incorporated  into  a  fully 
conservative,  time-dependent,  two-dimensional  code  in  order  to  test  its 
capability  in  predicting  subsonic  flows.  The  code  utilized  Van  Leer's^®  flux 
vector  splitting  technique  to  develop  a  finite  difference  scheme.  The  effort 
was  unsuccessful  in  that  the  model  did  not  cause  separation.  It  seemed  that 
the  smearing  inherent  in  the  scheme  was  the  cause  of  the  difficulty.  Unlike 
the  X  Scheme, the  conservative  scheme  spread  the  vortex  sheet  evenly  on  the 
surface  at  the  separation  point.  As  discussed  in  Section  2  of  this  report, 
the  X  Scheme  naturally  avoids  differentiating  the  entropy  and  in  some  sense 
the  velocity  at  a  point  where  the  velocity  changes  sign,  i.e.,  the  separation 
point.  This  property  seems  to  be  needed  in  order  to  force  separation  within 
the  model  used  here.  An  attempt  was  made  to  incorporate  the  model  into  a 
conservative  central  difference  code,  again  with  no  success.  Both  these 
schemes  are  notorious  for  spreading  convective  currents.  In  Ref  37  it  was 
found  that  both  these  schemes  were  inaccurate  for  boundary  layer  flows  because 
of  this  smearing. 

An  unsteady  x  Scheme  code  was  written  to  substantiate  the  suspicion 
stated  above.  The  separation  model  used  in  supersonic  flow  was  successful  in 
forcing  separation  on  a  cylinder  at  =  0.3  (nearly  incompressible).  The 
results  are  shown  in  Fig.  31.  Separation  was  forced  at  150“  from  the  forward 
stagnation  point,  as  can  be  seen  in  the  streamlines  of  Fig.  31a.  It  should  be 
noted  that  symmetry  was  assumed  in  the  computation.  The  surface  pressure  is 
shown  in  Fig.  31b.  These  results  indicate  that  the  model  can  be  used  in  the 
subsonic  region,  but  not  in  conjunction  with  a  conservative  formulation. 

The  original  intent  of  this  aspect  of  the  work  was  to  incorporate  the 
flow-separation  model  into  an  existing  subsonic,  conservative,  three- 
dimensional  code.  Since  the  available  codes  were  able  to  predict  flow 
separation  with  the  model,  in  its  present  form  this  aspect  of  the  work  was 
delayed  to  a  future  effort. 


8.  ASYMMETRIC  FLOWS 


I 

In  the  late  1940's,  the  phenomenon  of  lateral  instability  due  to 
asymmetric  vortex  formation  was  discovered.  On  slender-nosed  aircraft  flying 
at  a  high  angle  of  attack  (a  a  3  x  nose  half  angle),  the  flow  which  should  be 
I  symmetric  about  the  plane  containing  the  free  stream  velocity  vector  can  be 

very  asymmetric  due  to  the  vortices.  While  the  separation  points  are  slightly 
asymmetric,  the  resulting  vortices  are  very  asymmetric,  resulting  in  large 
side  forces  and  lateral  instability.  The  cause  of  this  asymmetry  has  been 
debated  since  it  was  first  discovered.  The  basic  issue  of  whether  it  is  a 
viscous  or  inviscid  phenomenon  is  still  uncertain.  The  inviscid  flow 
separation  model  used  here  is  an  ideal  tool  for  investigating  this  issue. 

One  of  the  first  indications  that  this  phenomenon  was  inviscid  in  nature 
came  from  the  work  of  Fiddes.^®  He  showed  small  disturbance  calculations  with 
vortex  sheets  modeled  by  vortex  filaments  that  exhibit  asymmetric  vortices. 

The  geometry  was  a  slender  cone  and  the  flow  was  incompressible.  Fiddes 
showed  that  two  solutions  existed  for  a  symmetrically  separated  flow.  One 
solution  being  symmetric,  the  resulting  vortices  were  symmetric,  and  the  other 
solution  asymmetric.  The  second,  asymmetric,  solution  exists  only  in  the  very 
high  angle-of-attack  regime,  a  >  2  x  cone  half  angle.  At  these  conditions, 
the  asytranetric  solution  compares  much  better  with  experimental  data  than  does 
the  symmetric  one.  Of  course,  in  the  wind  tunnel  experiment  Fiddes^®  found 
that  in  asymmetric  flow  the  separation  points  cannot  be  symmetric.  Very  small 
asymmetries  in  the  separation  locations  resulted  in  very  large  side  forces, 
indicating  very  large  asymmetries  in  the  vortex  locations.  Computationally, 
Fiddes  found  that  for  slightly  asymmetric  separation  point  locations  there 
existed  two  solutions,  one  which  was  slightly  asymmetric  and  the  other  very 
asymmetric.  Another  indication  that  the  phenomenon  is  inviscid  is  that 
experimental  data  indicate  an  insensitivity  of  the  onset  of  the  raise  in  side 
force  to  Reynolds  number. 

The  question  of  whether  the  small  disturbance  approximation  could  be  the 
cause  of  the  asymmetry  in  Fiddes  computation  can  be  answered  quite  simply  by 
the  Euler  computation  presented  here.  The  work  presented  here  is  confined  to 
the  supersonic  regime,  where  experimentalists  find  reduced  side  forces. The 
Euler  code  used  here  in  conjunction  with  the  flow-separation  model  was  able  to 


predict  asymmetric  flows.  Figure  32  shows  a  bifurcation  map  for  a  7”  cone  at 
=  1.6  with  separation  forces  symmetrically  at  9  =  25“  and  55“ .  The  plot 
shows  the  vertical  location  of  the  vortices  as  a  function  of  relative 
incidence.  Below  an  a/s  =  2.3,  only  one  solution  exists,  i.e.,  the  symmetric 
one.  When  the  angle  of  attack  is  increased  further,  two  solutions  exist:  the 
symmetric  and  the  asymmetric.  Note  the  three  vortex  locations  in  Fig.  32. 

The  two  asynmetric  vortex  locations  bracket  the  symmetric  one.  Unlike  Fiddes, 
finding  the  symmetric  solution  was  found  to  be  unstable  in  the  present 
computation.  The  only  way  to  achieve  a  symmetric  solution  was  to  enforce  a 
symmetry  condition.  If  the  symmetry  conditions  were  relaxed  after  a  symmetric 

solution  was  converged,  the  solution  could  revert  to  an  asymmetric  one. 

Figures  33  to  36  show  results  for  this  7”  cone  at  =  1.6  and 
a  =  20“.  Figure  33  shows  the  crossflow  streamlines  in  the  vicinity  of  the  lee 
plane;  the  asymmetry  of  the  vortices  is  obvious.  The  separation  points  are 

placed  symmetrically  (a  =  25“  and  155“)  as  is  the  case  in  all  the  calculations 

to  follow.  Figure  34  shows  the  surface  pressure  distribution;  it  should  be 
pointed  out  that  9  is  measured  counterclockwise.  The  crossflow  is  locally 
supersonic  so  that  shocks  precede  the  separation  points  (see  the  isobars  of 
Fig.  35).  The  vortex  to  the  left  is  much  closer  to  the  surface  (Fig.  33)  than 
that  on  the  right;  this  can  be  seen  in  its  influence  on  the  surface 
pressure.  The  minimum  surface  pressure  at  9  =  210“  is  a  result  of  this  vortex 
interaction.  The  isobars  of  Fig.  35  show  the  shocks  before  the  separation 
points  and  the  closed  isobars  representing  the  local  minimum  in  pressure  at 
the  vortex  centers.  While  the  asymmetry  is  significant,  it  is  localized  in 
supersonic  flow  because  of  the  limited  domain  of  influence.  Figure  34  shows 
that  the  surface  pressure  is  symmetric,  i.e.,  unaffected,  before  the 
separation  points.  This  is  the  reason  side  forces  due  to  asymmetries  decrease 
with  increasing  Mach  number.  The  asymmetry  may  be  significant,  but  because 
its  domain  of  influence  is  limited,  the  side  forces  may  be  small.  Figure  36 
shows  the  whole  crossflow  plane  with  the  streamlines;  note  that  the  flow  looks 
symmetric  except  in  the  vicinity  of  the  vortices. 

Figures  37  and  38  show  the  same  20“  angle-of-attack  case  with  symmetry 
enforced.  Figure  37  shows  the  crossflow  streamlines  and  Fig.  38  the  surface 
pressure  distribution.  Note  that  the  asymmetric  vortices  of  Fig.  33  bracket 
those  of  Fig.  37.  Figures  39  and  40  show  results  of  a  computation  at  small 
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enough  angles  of  attack  (a  =  15“ )  so  that  only  a  symmetrical  solution 
exists.  Figure  39  shows  the  crossflow  streamlines  and  Fig.  40  the  surface 
pressure.  It  should  be  pointed  out  that  this  calculation  was  performed  with 
no  symmetry  condition  imposed,  indicating  that  the  code  is  capable  of 
producing  symmetrical  results. 

Figure  32  indicates  that  the  flow  becomes  asymmetric  quickly  above 
a/«  =  2.3;  this  is  reinforced  in  Fig.  41.  The  figure  shows  the  crossflow 
streamlines  at  a  =  19“  (o/6  =  2.7)  to  be  very  asymmetric.  Figure  42  shows  the 
surface  pressure  which  seems  to  indicate  no  shock  at  separation  on  the  right 
side  and  a  small  shock  on  the  left  side.  Because  the  vortex  sheet  leaves  the 
body  at  such  a  large  angle  on  the  right  side  it  seems  the  flow  is  inhibited 
from  expanding  supersonically,  while  the  left-side,  shallow  vortex  allows  the 
flow  to  expand.  Figure  32  indicates  a  maximum  in  asymmetry  at  a  =  19“ . 

Figures  43,  44,  and  45  show  results  for  a  higher  angle  of  attack  a  =  23°, 
a/5  =  3.3.  While  the  vortices  seem  to  approach  the  same  height  above  the  cone 
(see  Fig.  43),  the  flow  is  by  no  means  approaching  symmetry.  This  is 
reinforced  by  the  surface  pressure  plot  of  Fig.  44.  It  seems  that  the  vortex 
height  may  not  be  the  best  indication  of  asymmetry.  The  asymmetries  in  this 
flow  is  further  evidenced  in  Fig.  45,  the  isobars. 

The  investigation  of  asymmetric  flows  is  incomplete  and  will  continue. 

The  ability  of  this  method  to  predict  the  onset  of  asymmetries  and  the  details 
of  these  flows  will  be  studied.  Computed  results  will  be  compared  with 
experimental  data.  The  method  presented  here  has  a  good  chance  of 
contributing  to  an  understanding  of  this  phenomenon.  It  seems  that  the 
interaction  of  the  two  vortices  is  the  cause  of  the  asymmetries;  this  will  be 
investigated. 
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9.  CONCLUSIONS 


The  main  conclusions  reached  as  a  result  of  this  work  are  listed  below: 

0  Shock-vorticity-induced  separation  and  shed-vorticity-induced 
separation  are  similar  in  that  the  vorticity  shed  into  the  flow 
field  is  reduced  smoothly  as  the  separation  point  is  moved  to  its 
shock-induced  location  (Fig.  9) 

0  Shock/vortex  interaction  can  cause  complex  shock  patterns  (Fig.  10 
and  11)  which  exhibit  interesting  fluid  physics 

0  The  relationship  between  sharp  leading-edge  separation  and  total 
pressure  losses  is  weaker  than  has  been  published  in  recent  years 

0  The  inclusion  of  a  proper  Kutta-type  condition  at  the  leading  edge 
of  a  sharp  leading-edge  wing  does  not  affect  results  very  much 
(Fig.  22) 

0  Secondary  separation  is  influenced  by  viscous  effects  more  than 
primary  separation  (Fig.  26) 

0  The  separation  model  was  simply  extendable  to  fully  three- 
dimensional  flow 

0  Difficulty  was  encountered  in  trying  to  apply  the  separation  model 
to  conservative  schemes 

0  Asymmetric  flow  can  be  predicted  with  the  separation  model  used 
here,  indicating  that  the  asymmetric  vortex  flow  phenomenon  is 
inviscid  in  nature. 
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Supersonic  Conical  Separation  Due  to  Shock  Vorticity 

Frank  Marconi* 

Grumman  Aerospace  Corporation,  Bethpage,  New  York 


A  numerical  study  of  Ihc  conical,  supersonic,  inviscid  flow  about  circular  slender  cones  at  large  incidence  is 
presented.  Although  the  flow  is  assumed  inviscid,  separation  and  an  accompanying  spiral  are  found.  This 
phenomenon  is  shown  to  be  a  valid  solution  to  Euler’s  epualions,  i.e„  valid  in  the  sense  of  being  independent  of 
any  numerical  viscosity.  The  separation  and  spiral  are  caused  solely  by  the  vorticity  produced  by  a  cross-flow 
shock.  The  resulting  separated  flowfields  arc  compared  to  full  potential  cakulations  and  to  espcrimcntal  data. 


Introduction 

H  REE-DIMENSIONAL  separated  flows  are,  in  general, 
accompanied  by  the  spiraling  of  streamlines  into  a  vortex. 
The  phenomenon  is  important  in  many  flow  situations:  the 
flow  about  the  leading  edge  of  a  wing;  the  flow  behind  bluff 
bodies;  and  the  flow  on  the  lee  side  of  slender  bodies  at  large 
angles  of  attack.  The  control  of  vortices  for  increased 
aerodynamic  efficiency  has  been  studied  by  aircraft  designers 
for  many  years  and  was  the  subject  of  a  recent  paper  by  Rao.' 
The  current  paper  is  a  numerical  study  of  spiraling  vortex 
flow  about  slender  cones  flying  at  supersonic  speeds  and  large 
incidence.  Although  the  cones  considered  here  may  be 
slender,  the  full  inviscid  Euler  equations  are  solved.  An 
understanding  of  these  flowfields  has  direct  applicability  to 
missile  aerodynamics  and  to  the  aerodynamics  of  aircraft 
forebodies. 

All  geometries  considered  in  this  numerical  study  are 
circular  cross-section  cones,  so  that  while  the  flow  Is  in  a  real 
sense  three-dimensional,  all  flow  variables  depend  on  only 
two  coordinates.  The  flowfield  is  assumed  supersonic 
everywhere  and  conical.  In  addition,  the  gas  is  assumed  ideal 
with  a  ratio  of  specific  heats  of  1.4  (i.e.,  air).  It  should  be 
realized  that  while  the  assumption  of  conical  flow  does 
simplify  the  numerical  problem,  the  conclusions  reached  have 
more  general  applicability.  The  basic  features  of  the  flowfield 
are  essentially  unchanged  in  the  fully  rhree-dimensional  flow. 

Conical  streamline  spirals  occur  in  the  separation  region  on 
the  leeward  side  of  slender  cones  at  large  angles  of  attack  (see 
Figs.  1  and  2).  The  existence  of  this  type  of  flow  pattern  in 
purely  inviscid  flow  has  been  the  subject  of  study  for  many 
years.  Kiichemann*  proved  analytically  that  for  two- 
dimensional  flow  with  distributed  vorticity,  inviscid 
separation  can  occur.  Two-dimensional  flow  separation  is 
accompanied  by  a  closed  recirculation  zone.  The  streamlines 
“trapped”  inside  this  recirculation  zone  cannot  be  tracked 
back  upstream  where  all  flowfield  conditions  are  known. 
Flow  properties  that  are  transmitted  along  streamlines  (for 
example,  entropy  and  vorticity)  must  somehow  be  determined 
inside  the  closed  recirculation  zone.  Fraenkel’  found  closed- 
form  solutions  for  the  incompressible  two-dimensional  flow 
over  cylinders  with  closed  recirculation  zones  caused  by 
distributed  vorticity  in  the  freestream.  Smith'  used  the 
analytical  work  of  Kiichemann  to  numerically  compute 
recirculation  zones  at  the  trailing  edge  of  a  two-dimensional 
airfoil.  The  vorticity  in  this  case  was  imposed  at  a  finite 
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distance  upstream  of  the  trailing  edge.  This  situation  is  similar 
to  vorticity  being  generated  by  a  transonic  shock  near  the 
trailing  edge  of  an  airfoil.  More  recently,  a  number  of 
researchers  working  on  solutions  to  Euler’s  equations  in  two- 
dimensional  transonic  flow’  *  have  found  inviscid  separation 
and  recirculation  zones  behind  cylinders  and  airfoils.  Salas’ 
showed  how  Euler  solutions  with  transonic  shocks  can  exhibit 
recirculation  zones  behind  cylinders.  The  vorticity  and  en¬ 
tropy  trapped  in  the  recirculation  zone  are  said  to  be  deter¬ 
mined  by  the  time-dependent  computation  used  to  asymptote 
to  a  steady  solution. 

In  fully  three-dimensional  flow  and  the  conical  flows 
considered  here,  separation  is  not  accompanied  by  a  closed 
recirculation  zone  but  by  spiraling  streamlines.’  This  point  is 
important  because  in  the  conical  flow  all  streamlines  originate 
in  the  freestream  where  all  flowfield  quantities  are  known. 
Therefore,  entropy  and  vorticity  can  be  uniquely  determined 
everywhere,  along  with  the  flowfield  and,  in  particular,  the 
required  shock  system.  Figure  1  shows  a  three-dimensional 
surface  streamline  starting  from  the  windward  plane  wetting 
the  body  surface  and  getting  swept  up  into  the  spiral.  Figure  2 
is  a  sketch  showing  how  the  cross-flow  streamlines,  which 
start  outside  the  bow  shock,  pass  through  the  shock  system, 
and  ultimately  get  swept  up  into  the  spiral  singularity  or  a  lee 
plane  node  (Fig.  2b).  The  dashed  line  of  Fig,  2b  divides  the 
flow  that  is  swept  into  the  spiral  node  from  that  which  goes 
into  the  lee  plane  node.  Cross- flow  streamlines  are  defined  as 
the  intersection  of  fully  three-dimensional  stream  surfaces 
and  the  surface  of  a  sphere  centered  at  the  apex  of  the  cone. 
Their  direction  in  the  x-y  plane  (Fig.  2)  can  be  defined  as  the 
ratio  of  the  two  components  of  the  cross-flow  velocity  u  and  f 
(Fig.  1). 

In  the  current  paper,  separation  will  be  loosely  defined  as 
ihe  process  by  which  a  stream  surface  originally  on  the  cone 
leaves  the  cone  surface  and  ultimately  forms  an  infinitely 
lurning  spiral  (see  Figs.  I  and  2).  Separation  occurs  along  a 
line  in  conical  flow;  ilie  intersection  of  this  line  and  the  cross- 
flow  plane  will  be  referred  to  as  the  separation  point.  In  a 
review  paper  on  conical  streamline  singularities,  Smiih' 
showed  analytically  that  a  conical  spiral  singularity  can  exist 
only  in  flows  with  vorticity.  In  particular,  he  found  that  a 
finite  radial  (spherical  r.  Fig.  1)  component  of  vorticity  was  a 
necessary  while  not  sufficient  condition  for  the  existence  of  a 
spiral  It  seems  clear  that  separation  and  an  accompaning 
spiral  are  linked  to  ihe  generation  of  vorticity  and  not  to  the 
details  of  the  viscous  processes  involved. 

There  are  three  ways  of  introducing  voriicity  in  an 
otherwise  irrotaiional  flow:  I)  by  specifying  vorticity  in  ihe 
freestream,  as  was  done  in  Refs.  3  and  4;  2)  through  shock 
waves;  and  3)  by  shedding  vorticily  from  ihe  boundary  layer 
ai  a  separation  point.  In  ihis  paper  only  the  vorticity  produced 
by  the  shock  system  will  be  discussed  m  detail.  It  should  be 
noted  that  currently  only  through  solutions  of  the  Euler 
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equations  can  this  aspect  of  the  flow  be  modeled  (i.e.,  full 
potential  shocks  do  not  generate  vorticity).  The  third 
mechanism  of  vorticity  introduction,  shedding,  has  been 
studied  extensively  by  Smith*'"’  and  Fiddes."  The  current 
author  has  shown  preliminary  Euler  solutions  that  use  this 
model  to  force  separation  at  a  specified  location. 

In  order  to  shed  vorticity,  the  flow  is  modeled  by  assuming 
that  there  exists  an  infinitesimally  thin  sheet  of  infinite 
vorticity  across  which  there  is  a  jump  in  velocity.  The  rest  of 
the  flow  can  be  assumed  irrotational.  The  vortex  sheet  will 
generally  roll  up  into  a  spiral.  The  sheet  is  forced  to  originate 
at  a  separation  point  by  imposing  a  Kutta  condition.  At  sharp 
edges  (sharp  wing  tip,  for  example),  the  condition  is  similar  to 
that  impost  at  wing  trailing  edges  in  three-dimensional  flow 
to  exclude  singularities  due  to  flow  turning  around  sharp 
edges.  For  smooth  bodies,  the  proper  Kutta  condition  to  be 
imposed  to  start  a  vortex  sheet  in  an  otherwise  completely 
irrotational  flow  has  been  studied  in  detail  in  Refs.  9  and  10. 
The  location  of  the  origin  of  shedding  on  the  body  is  the  only 
free  parameter  and  is  specified  to  match  the  boundary-layer 
separation  point.  In  contrast  to  bodies  with  salient  edges, 
where  the  edge  is  taken  as  the  separation  point,  smooth  bodies 
require  that  the  separation  point  be  computed  from  a 
boundary  layer  coupled  to  an  inviscid  flow  iteration"  or  from 
experimental  data.  Calculations  including  forced  separation 
from  a  shedding  vortex  sheet  have  been  performed  for  a 
number  of  years  (for  example,  the  early  work  of  Bryson"). 
More  recently,  Rizzi"  has  shown  solutions  to  Euler’s 
equations  with  separation  from  sharp  leading  edge  wings,  and 
Nielsen"  has  shown  solutions  to  Euler's  equations  with 
separation  from  smooth  bodies.  It  should  be  pointed  out  that 
the  work  of  Refs.  14  and  IS  seem  inconsistent  with  the 
analysis  of  Refs.  9  and  10  for  the  flow  near  the  origin  of  this 
sheet  (i.e.,  the  separation  point),  while  the  present  author  has 
used  the  model  of  Ref.  9  to  force  separation  at  a  specified 
location." 

It  should  be  made  clear  at  this  point  that  the  two 
mechanisms  of  vorticity  introduction,  shock  waves  and 
shedding,  both  produce  valid  solutions  to  Euler’s  equations 
with  separation  and  spirals.  While  the  two  mechanisms  pro¬ 
duce  similar  flowfield  structures,  the  details  of  the  separation 
and  spiraling  sheet  are  different.  The  differences  in  detail  are 
outlined  in  Ref.  12.  The  main  difference  should  be  clear.  The 
separation  caused  by  the  shock  vorticity  occurs  with  no  ad¬ 
ditional  conditions  and  uniquely  determines  separation  point 
location,  while  the  forced  separation  point  location  is  ar¬ 
bitrary.  It  is  the  arbitrariness  of  the  separation  point  location 
that  allows  the  freedom  necessary  to  match  the  boundary- 
layer  separation  point  location.  The  separation  point  in  the 
shock-induced  case  may  bear  little  relation  to  the  real  viscous 
flow,  but  the  shock-produced  vonicity  can  play  a  role  in  the 
structure  of  the  real  separated  flow.  The  point  here  is  that  the 
vorticity  in  the  flowfield  comes  from  two  sources — shock 
waves  and  shed  from  the  boundary  layer — both  of  which 
contribute  to  the  quantitative  structure  of  the  flow.  Only  the 
shock  vorticity  will  be  considered  herein,  so  it  is  not  surprising 
that  the  comparisons  with  experimental  data  are  poor  in  some 
sense;  in  particular,  the  separation  point  locations  are 
significantly  different. 

The  intent  of  this  paper  is  to  show  that  the  vorticity 
produced  by  shock  waves  can  cause  separation  and  that 
solutions  to  the  inviscid  Euler  equations  exhibiting  this 
phenomenon  are  valid,  in  the  sense  of  being  independent  of 
any  numerical  artifices.  The  qualitative  structure  of  the 
separated  flow  produced  by  shock  vorticity  will  be  shown  to 
be  the  same  as  that  produced  in  the  presence  of  a  boundary 
layer,  indicating  the  possibility  of  this  vorticity  being 
significant,  in  the  section  that  follows,  the  overall  com¬ 
putational  procedure  will  be  outlined.  The  issues  of  con¬ 
vergence  both  in  terms  of  iteration  scheme  and  grid 
refinement  will  be  addressed  in  this  section.  Then  a  detailed 
study  of  computational  results  will  be  presented  including 


comparisons  with  experimental  data  and  full  potential 
calculations.  The  final  section  of  the  paper  will  present  a 
summary  of  the  findings  of  this  work  and  the  future  direc¬ 
tions  of  this  effort. 

Computational  Procedure 

The  overall  numerical  procedure  used  in  this  study, 
although  employing  a  number  of  new  features,  is  essentially 
standard."  The  fully  three-dimensional  Euler’s  equations  are 
solved  with  an  explicit  marching  technique.  The  marching 
direction  (z.  Fig.  I)  is  an  iterative  coordinate  for  the  conical 
flow  considered  here.  The  scheme  is  restricted  by  the  fact  that 
the  axial  component  of  velocity  (fV,  Fig.  I)  must  be  super¬ 
sonic  everywhere.  The  marching  scheme  is  continued  until  the 
flowfield  is  invariant  with  respect  to  the  computational 
marching  direction  except  for  a  scale  factor.  The  finite- 
difference  scheme  used  is  Moretti’s  characteristic  based  X- 
scheme. All  shock  waves  in  these  flowfields  are  fit  and  are 
forced  to  satisfy  the  exact  Rankine-Hugoniot  jump  con¬ 
ditions.  The  bow  shock  is  fit  as  the  outermost  boundary  of  the 
flowfield.  On  the  low  pressure  side  of  the  bow  shock, 
freestream  conditions  exist.  The  cross-flow  shock  is  fit  as  an 
internal  boundary  of  the  flowfield  with  its  low  pressure  side 
being  computed  as  the  computation  proceeds. " 

The  cross-flow  shock  computation  is  a  critical  part  of  the 
overall  procedure,  particularly  since  it  plays  such  a  critical 
role  in  separation.  As  mentioned  in  the  previous  section,  it  is 
the  vorticity  produced  by  the  cross-flow  shock  that  produces 
the  spiral.  Any  scheme  that  captures  the  cross-flow  shock  and 
introduces  additional  artificial  viscosity  to  stabilize  it  runs  the 
risk  of  distorting  the  separation  and  spiral.  The  bow  shock 
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Fig.  1  ThreeHliineiisioiitl  sketch  of  flowfieM  and  coordinate  svstem. 
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also  iniroduces  vorticity  in  the  nowHeld,  but  not  enough  to 
produce  separation.  The  cross-flow  shock  is  caused  by  the 
fact  that  at  a  large  incidence  the  9  component  of  the  velocity 
(u.  Fig.  I)  becomes  supercritical.  The  flow  expands  in  going 
from  the  windward  symmetry  plane,  and  u  can  become 
supercritical  if  the  incidence  is  high  enough.  This  component 
of  velocity  must  vanish  in  the  lee  symmetry  plane  due  to  the 
boundary  conditions.  A  supercritical  cross  flow  generally 
passes  through  a  relatively  normal  shock  before  stagnating  in 
the  lee  plane.  The  cross-flow  shock  exists  to  turn  the  three- 
dimensional  flow  parallel  to  the  lee  symmetry  plane.  The 
variation  in  strength  of  the  cross-flow  shock  can  be  quite 
large.  The  shock  is  strongest  at  the  body  where  u  is  highly 
supersonic  and  approaches  zero  strength  in  the  field  where  u 
becomes  sonic.  It  is  this  variation  in  strength  that  produces  a 
cross-flow  entropy  gradient  and  thus  radial  vorticity.  This 
vorticity.  in  turn,  causes  shock-induced  separation. 

The  low-pressure  side  of  the  cross-How  shock  is  computed 
with  one-sided  differences  away  from  the  shock.  This  is 
consistent  with  the  fact  that  the  cross  flow  is  supersonic.  The 
Rankine-Hugoniot  conditions  and  the  compatibility  condition 
along  a  bicharacteristic  reaching  the  shock  on  its  high- 
pressure  side  supply  enough  information  to  compute  the 
deviation  of  each  shock  point  from  a  conical  ray  in  addition 
to  all  the  primitive  variables  on  its  high-pressure  side.  The 
bicharacteristic  used  is  the  one  in  the  plane  containing  the 
local  normal  to  the  shock  and  the  marching  direction.  All 
shock  points  are  computed  with  the  postcorrection  scheme 
proposed  by  Rudman'’  and  independently  by  deNeef.“  The 
shocks  converge  (i.e.,  each  shock  point  becomes  aligned  with 
a  conical  ray)  with  the  rest  of  the  flowfield.  In  all  the  com¬ 
putations  presented  here,  the  last  cross-flow  shock  point  fit 
had  a  normal  Mach  number  of  about  1.05  (pressure  ratio  of 
1.12).  The  finite  difference  scheme  was  able  to  capture  weaker 
shock  points. 

.An  exponential  stretching  is  used  to  cluster  grid  points  near 
the  surface  of  the  cone.  This  is  important  because  flow 
conditions  that  exhibit  the  spiral  singularity  also  exhibit  large 
gradients  near  the  body  behind  the  cross-flow  shock.  Without 
the  stretching,  the  radial  resolution  behind  the  shock  can  be 
too  coarse  to  pick  up  the  separation  and  spiral  (see,  for 
example,  the  Euler  calculations  made  using  a  procedure  very 
similar  to  the  one  used  here  without  a  stretching  present  in 
Ref.  21).  The  stretching  is  of  the  form: 

r=  -t„(/~(/-e-'^)X)(C-B)/iS  +  B  (I) 

where  r  is  the  polar  radius  in  thex-y  plane  (not  to  be  confused 
with  the  spherical  radius  r  of  Fig.  I),  and  X  is  the  com- 
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putational  coordinate  going  from  0  on  the  body  to  I  at  the 
bow  shock.  The  other  computational  coordinate  is  V  .  which  is 
a  linearly  normalized  polar  angle  9.  Y  is  normalized  to  go 
from  0  in  the  windward  symmetry  plane  to  1  on  the  low- 
pressure  side  of  the  cross-flow  shock  and  from  0  again  on  the 
high-pressure  side  of  the  shock  to  1  in  the  lee  plane.  The 
transformation  to  the  final  computational  plane  is  completed 
by  Z=;.  In  the  transformation.  Hq.  (I).  Sand  Care  the  polar 
radii  of  the  body  and  bow  shock,  and  d  is  a  clustering 
parameter.  The  grid  points  are  evenly  spaced  in  each  .V-  > 
plane.  The  transformation,  Eq.  (I),  clusters  grid  points  near 
the  body  depending  on  the  value  of  d  used.  The  trans¬ 
formation  approaches  linear  as  d  goes  to  zero.  All  of  ihe 
calculations  shown  here  were  performed  with  d  =  2.5.  The 
gradient  in  the  strength  of  the  cross-flow  shock  is  quite  large, 
which  is  the  reason  for  the  large  vorticity  produced  behind  it. 
The  clustering  of  grid  points  near  the  body  allows  for  enough 
resolution  to  compute  this  shock  decay  properly.  The  number 
of  points  required  to  resolve  this  phenomenon  without  a 
stretching  would  be  inordinately  large. 

Shock-induced  inviscid  separation  produces  a  contact 
surface  emanating  from  the  separation  point  (Fig.  2).  It  is  this 
contact  that  ultimately  spirals.  The  contact  surface  has  a  jump 
in  entropy  and  an  accompanying  jump  in  velocity.  In  conical 
flow,  entropy  is  constant  on  cross-flow  streamlines.  With  this 
fact,  one  can  see  (Fig.  2)  that  the  entropy  on  the  wind  side  of 
the  contact  comes  from  the  wind  symmetry  plane  bow  shock. 
The  streamline  wets  the  body  and  passes  through  the  base  of 
the  cross-flow  shock.  This  is  the  high  entropy  side  of  the 
contact.  There  are  two  possibilities  for  the  entropy  on  the  lee 
side  of  the  contact.  The  one  shown  in  Fig.  2a  corresponds  to 
the  lee  bow  shock  entropy  wetting  the  lee  side  of  the  cone  and 
ending  up  on  the  lee  side  of  the  contact.  There  is  a  possibility 
of  a  saddle  point  existing  in  the  lee  symmetry  plane  off  the 
body  as  shown  in  Fig.  2b.  In  this  case,  the  entropy  on  the  lee 
side  of  the  contact  would  come  from  this  saddle.  In  either  case 
the  cross-flow  steamlines  that  wet  the  body  are  tracked  in  each 
step  of  the  iteration,  and  the  proper  entropy  is  imposed  on  the 
cone  including  the  jump  in  entropy  at  the  separation  point. 
The  entropy  discontinuity  off  the  body  is  captured  with  the 
finite  difference  scheme. 

The  finite  difference  scheme'’  used  requires  that  all  spatial 
differences  (computational  X  and  Y)  be  directional, 
depending  on  characteristic  directions.  In  the  case  of  entropy, 
the  characteristic  direction  is  the  particle  path  or  streamlines, 
so  that  both  X  and  Y  derivatives  of  entropy  are  taken 
depending  on  the  direction  of  the  local  flow.  The  scheme 
basically  avoids  differences  across  the  contact.  In  order  to 
achieve  formal,  second-order  accuracy,  the  X-scheme  uses 
three-point,  one-sided  differences.  In  this  flowfield,  the  three- 
point  differences  in  entropy  sometimes  required  crossing  the 
contact,  which  caused  wiggles.  In  order  to  ensure 
monotonicity  at  the  sheet,  these  differences  in  entropy  were 
taken  with  two-point  windward  formulas.  The  sheet  was 
captured  between  two  points  until  it  began  to  roll  up.  The 
region  where  the  sheet  rolls  up  is  somewhat  smeared,  but  the 
scheme  seems  to  describe  this  region  adequately.  In  any  in¬ 
viscid  calculation,  approximations  must  be  made  for  the 
region  near  the  center  of  the  spiral.  In  the  irrotational 
calculations  of  Refs.  II  and  13,  the  center  of  the  spiral  is 
replaced  by  a  concentrated  vortex.  One  advantage  of  solving 
the  Euler  equations  is  that  this  region  of  concentrated  vor- 
ticiiy  can  be  computed  automatically 

It  was  proven  by  Smith^  that  a  forced  sheet  leaves  a  smooth 
body  tangent  to  the  surface  and  that  there  must  be  a  jump  in 
cross-flow  velocity  at  separation.  The  cross-flow  velocity  on 
the  windward  side  of  the  sheet  is  finite,  while  the  cross  How  on 
the  lee  side  stagnates.  The  arguments  of  Ref  f  are  based  on 
irrotational  and  iseniropic  assumptions  and  do  not  hold  lor 
ilows  with  cross-flow  shocks.  It  can  be  shown  that  when  the 
sheet  exhibits  a  jump  in  emropy.  the  cross-llow  velocines  on 
cither  side  of  the  sheet  can  match  and  still  mainiain  con 
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tinuous  pressure  across  the  sheet.  Of  course  the  modulus  of 
the  velocity  must  then  jump.  This  jump  is  absorbed  in  the 
radial  component  of  velocity  (w.  Fig.  1).  This  seems  to  be  the 
situation  in  the  inviscid  shock-induced  separation.  Results 
will  be  shown  later  that  indicate  the  continuous  variation  of 
cross-flow  through  the  separation  point  for  the  shock-induced 
case. 

In  order  to  speed  convergence,  grid  refinement  was  used  in 
all  the  calculations  shown  in  this  paper.  The  calculations  were 
started  with  a  10x10  (polar  radial  x  circumferential)  grid. 
The  calculation  was  continued  until  a  reasonable  convergence 
was  reached,  and  then  the  grid  was  halved.  The  grid  was 
refined  three  limes,  the  final  grid  having  73  x  73  points  in  each 
X-  Y  plane.  The  iteration  was  continued  on  the  finest  grid  until 
the  maximum  variation  in  pressure  in  the  spherical  radial  (r. 
Fig.  1)  direction  was  less  than  1  x  10~\  This  is  the  maximum 
correction  in  pressure.  The  solutions  are  essentially  converged 
at  a  maximum  correaion  of  10  so  that  10  is  a  con¬ 
servative  convergence  criterion.  This  is  a  reduction  of  about 
nine  orders  of  magnitude  from  the  beginning  of  the  iteration. 
Critical  parameters  such  as  shock  location  and  stagnation 
point  location  are  constant  to  five  places  by  the  time  the 
calculation  is  stopped. 

The  sensitivity  of  this  calculation  to  the  initial  conditions 
used  was  studied  early  in  this  work.  It  was  found  that  the 
captured  shock,  coarse  grid,  initial  conditions  to  the  fine  grid 
exhibited  a  spiral.  The  spiral  was  eliminated  before  starting 
the  shock  fit,  fine  grid  calculations  with  no  change  in  the  final 
result.  The  spiral  simply  developed  again.  In  a  numerical 
experiment,  a  converged  high  a  (with  spiral)  result  was  used  as 
initial  condition  for  an  a  =  0  case.  The  incidence  was  reduced 
gradually,  and  the  cross-flow  shock  was  weakened  and 
ultimately  was  eliminated  along  with  the  spiral.  The  opposite 
sequence  was  tried  and  the  spiral  reappeared.  The  numerical 
experiments  Indicated  an  independence  of  the  converged 
solution  to  initial  conditions. 

The  finite  difference  scheme  used  in  this  work  is  an  explicit 
marching  scheme  that  is  notoriously  inefficient  for  converging 
to  a  conical  or  steady  solution.  A  number  of  speed-up 
procedures  were  tried,  including  multigrid  and  local  marching 
steps.’  None  of  these  worked.  Local  marching  steps  did 
produce  a  significant  reduction  in  running  times,  but  proved 
unstable  in  the  presence  of  the  cross-flow  shock.  The  cross- 
flow  shock  was  captured  with  the  X-scheme  in  the  first  two 
coarse  grid  calculations.  The  shock  was  detected  and  fit  from 
the  third  grid  on.  While  the  captured  shock  location  was 
significantly  off  the  converged  result,  this  scheme  eliminated 
some  shock  detection  problems.  One  advantage  of  the  explicit 
marching  scheme  used  was  that  It  is  totally  "vectorizable." 
The  time  consuming  parts  of  the  code  (i.e.,  interior  point 
computation)  utilized  the  vector  architeaure  of  the  CRAY  I 
computer.  This  made  the  computation  about  20  times  faster 
on  the  CRAY  than  on  an  IBM  3033.  All  the  computations 
shown  in  this  report  took  about  30  CPU  minutes  on  the 
CRAY  1. 

Figure  3  shows  a  comparison  of  the  surface  pressure  for  the 
flow  field  about  a  10  deg  half-angle  cone  at  =  2  and  a  =  25 
deg.  The  symbols  are  the  pressure  computed  with  the  standard 
73  X  73  grid.  The  solid  line  was  computed  with  a  finer  89  x  89 
grid.  The  figure  shows  that  the  fine  and  standard  grid  surface 
pressures  are  virtually  identical,  indicating  that  the  standard 
73  x  73  grid  resolves  all  flowfield  features.  In  both 
calculations,  about  35%  of  the  circumferential  points  are 
behind  the  cross-flow  shock  (i.e..  27  points  for  the  73x73 
grid).  This  flowfield  has  a  shock-induced  separation  at 
149  6  deg,  indicated  by  the  minimum  in  surface  pressure 
behind  the  shock  ($•  165  deg).  This  minimum  in  pressure  is 
just  below  the  core  of  the  spiral.  Both  computations.  7)  x  73 
and  89  x  89.  give  the  same  separation  point  and  spiral 
location,  indicating  an  independence  of  ihe  results  to 
numerical  viscosity  due  to  truncation  error  As  a  matter  of 
fact,  ihis  flowfield  was  computed  with  an  even  coarser  grid 


(37  X  37)  that  exhibits  the  same  basic  features  (i.e..  separation 
and  spiral).  While  it  is  very  difficult  to  prove  the  existence  of 
shock-induced  inviscid  separation,  the  aforementioned 
numerical  experiments  Indicate  the  validity  of  such  solutions. 
Additionally,  the  analytical  work  for  two-dimensional  flow 
reported  in  Refs.  2  and  3  confirm  the  possibility  of  such 
solutions.  The  uniqueness  of  the  conical  flow  solution  would 
seem  to  be  on  firmer  ground  because  of  the  absence  of  closed 
recirculation  zones. 

Computational  Results 

In  this  section,  computational  results  showing  spiraling 
streamlines  due  to  the  vorticity  produced  by  a  cross-flow 
shock  will  be  analyzed.  The  first  configuration  considered  is  a 
10-deg  cone  at  =2  and  a  =  25  deg.  The  surface  pressure 
distribution  has  been  shown  in  Fig.  3.  The  cross-flow 
streamlines  are  shown  in  Fig.  4,  along  with  the  computed  bow 
and  cross-flow  shocks.  The  streamlines  show  saddles  in  the 
wind  and  lee  symmetry  planes  on  the  cone  surface.  In  ad- 
dition,  there  is  a  saddle  at  the  separation  point  and  a  node  at 
the  center  of  the  spiral  (this  configuration  is  sketched  in  Fig. 
2a).  The  separation  point  and  spiral  can  be  seen  more  clearly 
in  the  blowup  of  Fig.  5.  If  the  nodes  and  saddles  are  summed 
using  the  procedure  of  Ref.  8,  it  can  be  shown  that  the  proper 
number  of  nodes  and  saddles  exist  in  the  streamline  pattern  of 
Fig.  4.  The  separation  point  (S,  =  149.6  deg)  corresponds  to 
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Fig.  6  SurfKt  croM-flow  vdocily  (M.  >  2,  i  =  10  deg,  a  -  25  deg). 


Fig.  ^  Isoban  in  cross  plane 
»  2.  &-  10  deg.  an  25  deg). 


the  plateau  in  pressure  (Fig.  3)  just  after  the  shock.  If  one 
considers  the  momentum  equation  in  the  D  direction,  it  in¬ 
dicates  that  dp/ 96  =  0  at  a  cross-flow  stagnation  point 
(u  =  c  =  0).  The  separation  point  in  this  flowfield  is  a  real 
cross-flow  stagnation  point.  The  streamlines  (Fig.  5)  show 
that  the  How  moves  in  the  negative  6  direction  from  the  lee 
stagnation  point  toward  the  separation  point.  In  doing  so.  the 
flow  expands  (i.e.,  drop  in  pressure  between  8=  180  deg  and 
165  deg).  The  flow  then  recompresses  to  the  separation 
point.  This  recompression  phenomenon  is  the  cause  of 
secondary  separation.  The  inviscid  separation  point  location 
IS  tar  from  that  of  the  viscous  flow.  This  can  be  surmised  by 
I  he  fact  that  the  shock  is  strong  enough  to  separate  a  bound¬ 
ary  layer  at  ns  base.  The  inviscid  separation  point  is  too  far 
downstream  of  the  shock.  The  streamlines  of  Fig.  5  clearly 
show  how  all  of  the  flow  is  ultimately  swept  up  into  the  in¬ 
finitely  turning  spiral.  The  apparent  power  of  the  Euler 
equations  to  describe  the  region  near  the  center  of  the  spiral 
should  be  noted.  A  detailed  study  of  the  streamline  pattern 
near  the  center  of  the  spiral  revealed  that  the  streamlines 
asymptote  to  an  ellipse.  This  was  first  indicated  by  Smith.'*  It 
also  should  be  pointed  out  that  Fig.  5  shows  that  streamlines 
wrapping  around  the  spiral  very  rapidly  approach  the  lee  side 
1)1  the  separation  line  (contact)  Consider,  for  example,  the 
ihird  sireamline  off  the  cone  in  Fig  5  It  wraps  around  ihe  top 
ol  the  cpiral.  comes  back,  and  approaches  the  contact  a  short 
distance  olf  the  bodv,  creating  an  cniropy  layer  on  the  lee  side 
111  the  coniaci.  This  eniropy  layer  lends  lo  weaken  ihc  coniaci 


Fig.  g  Lines  of  constant  radial  vorticiiy  x  spherical  radius  i  Vf,  =  2. 
5  ^  10  deg.  a  >  25  deg). 


Figure  5  also  shows  that  the  separation  line  leaves  the  cone 
at  a  finite  angle.  Smith's  analysis^  concluded  that  at  a  forced 
separation  the  sheet  comes  off  tangent.  However,  this  dif¬ 
ference  is  due  to  the  fact  that  the  separation  in  the  shock- 
induced  case  is  nonisentropic  and  is  a  real  cross-flow 
stagnation  point.  Figure  6  shows  the  cross-flow  velocity  on 
the  surface  of  the  body  (u  =  0  from  the  body  boundary 
condition).  The  point  to  be  considered  here  is  that  the  velocity 
passes  through  zero  (separation)  smoothly.  The  analysis  of 
Smith  indicates  that  separation  occurs  at  a  discontinuity  in 
cross  flow,  the  jump  in  velocity  determining  the  sheet 
strength.  This  is  a  basic  difference  between  forced  and  shock- 
induced  separation  and  is  discussed  in  detail  in  Ref.  12.  The 
jump  in  cross-flow  velocity  in  Fig.  6  is  due  to  the  shtxk.  It  is 
also  interesting  to  note  the  maximum  negative  u  that  occurs 
under  the  spiral.  This  negative  cross  flow  can  become 
supersonic,  causing  a  second  reverse  cross-flow  shock.  A 
sample  of  this  will  be  shown  later. 

Figure?  shows  the  isobars  for  the  same  case  (M,  =  2,  6=  10 
deg,  a  =  25  deg).  It  should  be  noted  how  smoothly  the  cross¬ 
flow  shock  transitions  to  zero  strength  in  the  field.  The  most 
interesting  aspect  of  the  figure  is  the  closed  isobar  at  the 
center  of  the  spiral.  It  represents  an  absolute  minimum  in  the 
flowfield  pressure.  The  component  of  vorticiiy  In  the 
spherical  radial  direction  is  given  by: 

/  du  dv  .  I 

W,  =  -  I  ucos*  ■)-  —  sino  -  —  )  /  ( rsino) 

'  d0  d6  '/ 

where  u  and  v  are  the  cross-flow  velocities  defined  in  Fig.  1 
and  r.  9,  0  are  spherical  coordinates  (Fig.  I).  Figure  8  shows 
lines  of  constant  K),.  The  figure  shows  how  the  vorticiiy  is 
distributed  and  is  produced  by  the  cross-flow  shock.  Note  that 
{),  IS  small  until  the  cross-flow  shock,  indicating  that  the  bow 
shock  doesn’t  produce  enough  vorticiiy  to  cause  separation  on 
Its  own.  The  vorticiiy  produced  by  the  bow  shock  is  an  order 
of  magnitude  smaller  than  that  produced  by  the  cross-flow 
shock  and  two  orders  of  magnitude  smaller  than  that  at  the 
center  of  the  spiral.  For  this  reason,  there  are  no  contours  of 
vorticiiy  outside  the  region  behind  the  cross-How  shock  m 
Fig  8.  The  vorticity  is  negative  causing  the  counterclockwise 
spiraling  of  this  flow,  and  us  absolute  magnitude  is  maximum 
at  the  center  of  the  vortex. 

The  results  shown  in  Figs  ')-ll  are  all  tor  the  same  flow 
conditions  (  Vf,  =  2.  6-  10,  i  =  20  Jegl  This  is  a  lower  angle 
ol  attack  than  the  previous  case,  so  that  the  spiral  is  smaller 
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and  has  moved  toward  the  lee  plane.  This  can  be  seen  by 
comparing  Figs.  5  and  10.  It  also  can  be  seen  from  the 
streamlines  (Figs.  9  and  10)  that  this  configuration  has  an 
additional  saddle/node  combination  in  the  lee  plane  as 
sketched  in  Fig.  2b.  The  dashed  line  ending  at  the  saddle  in  the 
lee  plane  divides  the  flow,  which  is  ultimately  swept  into  the 
spiral  node  from  that  which  ends  up  in  the  node  in  the  lee 
symmetry  plane.  The  entropy  on  the  dashed  streamline 
ultimately  wets  the  lee  side  of  the  contact  (separation  line). 
This  typte  of  structure  has  been  noticed  experimentally."  The 
phenomenon  is  similar  to  the  standard  “liftofP*  conical  flow 
structure  noticed  at  lower  incidence  with  no  spiral.  In  the 
present  study,  subcritical  incidences  were  computed  for  the 
same  lO-deg  cone  at  =  2,  and  liftoff  was  observed  with  no 
vortex.  As  the  incidence  is  raised,  the  flow  becomes  super¬ 
critical  and  the  spiral  begins  to  form  near  the  lee  plane.  The 
liftoff  remains,  as  shown  in  Figs.  9  and  10.  As  the  incidence  is 
raised  still  further,  the  spiral  gets  larger  and  the  liftoff  is 
eliminated  (see  Figs.  4  and  3).  It  should  be  noted  that  the 
limiting  streamline  (dashed  line  in  Fig.  10)  seems  not  to  ap¬ 
proach  the  saddle  normally.  This  is  possible  if  distributed 
vorticity  is  present  in  the  flow,-^  which  is  the  case  here.  The 
isobars  in  Fig.  1 1  show  that  the  vortex  is  so  weak  that  the 
minimum  in  pressure  at  its  center  does  not  appear  due  to  the 
resolution  of  the  contours. 

Figure  12  is  a  plot  of  the  inviscid  separation  point  location 
vs  a  measure  of  shock  strength.  The  cone  has  a  lO-deg  half- 
angle  and  the  freestream  Mach  number  is  2.  The  cross-flow 
shock  strength  was  increased  by  increasing  incidence.  The 
pressure  ratio  used  is  the  one  at  the  base  of  the  shock;  of 
course,  all  cross-flow  shocks  go  to  zero  strength  at  their  tips. 
The  pressure  ratio  used  is  a  good  measure  of  the  shock 
strength  variation  and  thus  the  vorticity  produced  by  the 
shock.  At  each  data  point,  the  corresponding  a  is  noted.  The 
plot  shows  how  as  the  shock  strength  decreases,  the 
separation  point  moves  to  the  lee  symmetry  plane.  At  the 
same  time,  the  spiral  region  is  getting  smaller.  The  highest 
incidence  computed  was  25  deg,  since  above  that  value  the 
axial  Mach  number  on  the  cone  surface  approaches  sonic, 
making  it  impossible  to  march.  At  the  low  end,  a  =  19  deg,  the 
spiral  was  so  small  that  it  was  difficult  to  resolve  numerically, 
and  no  lower  incidences  are  shown.  The  interesting  feature  of 
this  figure  is  that  an  extrapolation  of  the  curve  would  seem  to 
indicate  that  inviscid  separation  moves  to  the  lee  plane  before 
the  shock  is  eliminated  (Pi/Pi  =  I.).  An  extrapolation  would 
indicate  that  the  spiral  is  eliminated  at  al.SS.  This 

corresponds  to  a  normal  Mach  number  slightly  above  1.3. 
which  is  approximately  the  region  where  the  full  potential 
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approximation  is  valid  for  these  flowfields.  It  would  seem 
that  below  a  maximum  normal  Mach  number  of  1.3,  the 
cross-flow  shock  may  not  produce  enough  vorticity  to  cause 
separation. 

Figures  13-15  deal  with  an  interesting  case(Af_  =  3,  4  =  9.46 
deg,  and  a  =  25  deg).  Figure  13  show  the  streamlines  and  the 
crcss-How  shock  near  the  lee  plane.  The  shock  exhibits  a  kink 
as  it  passes  from  the  influence  of  the  spiral.  Near  the  cone,  the 
shock  must  deflect  the  streamlines  upward  in  order  that  they 
may  pass  over  the  spiral;  beyond  the  top  of  the  spiral,  this  is 
no  longer  true  and  the  shock  acts  like  a  normal  shock.  The 
two  regions  are  separated  by  the  kink  in  the  shock,  and  the 
shock  slope  in  the  cross-flow  plane  had  to  be  differenced  away 
from  this  point.  Another  interesting  feature  of  this  flow  is  the 
fact  that  the  expansion  of  the  reverse  flow  from  the  lee 
stagnation  point  is  so  large  that  the  negative  cross  flow 
becomes  supersonic  near  the  body.  The  smooth  recom¬ 
pression  shown  in  Fig.  3  is  replaced  by  a  reverse  cross-flow 
shock.  This  phenomenon  has  been  noted  experimentally  in 
Ref.  24.  The  reverse  cross-flow  shock,  which  was  captured, 
can  be  seen  in  the  isobars  of  Fig.  14.  The  second  (reverse) 
shock  is  on  the  lee  side  of  the  primary  cross-flow  shock.  It  is 
indicated  by  the  clustering  of  the  isobars  between 
i'n(p/p,  )  =  -  1.2  and  f«(p/p,)=  -0.73  on  the  cone  surface. 
The  shock  is  not  strong  enough  to  produce  a  secondary  in- 


Fig.  13  Crou-flow  simmliacs  near  lee  plane  =  3,  i  =  9.46  deg, 
a  =  IS  deg). 


Mg.  14  Isoban  near  lee  plane  (Af,  «  3.  i  ■  9.46  deg,  a  >  25  deg). 


viscid  separation,  whereas  in  the  experiment  of  Ref.  24  it  was 
strong  enough  to  separate  the  boundary  layer.  The  possibility 
of  having  secondary  separation  produced  solely  by  shock 
vorticity  is  an  interesting  topic  for  future  study.  The  strength 
of  the  vortex  is  also  indicated  in  Fig.  14  by  the  closed  isobars 
representing  a  steep  pressure  minimum  at  the  vortex  center. 

Figure  15  shows  the  surface  pressure  for  the  .Vf,  =3, 
4  =  9.46  deg,  and  a  =  25  deg  case.  The  solid  line  is  the  present 
calculation  with  shock  vorticity-induced  separation  at 
4=  155.3  deg,  again  at  the  plateau  in  pressure  just  after  the 
primary  cross-flow  shock.  The  reverse  cross-flow  shock  can 
be  seen  in  the  pressure  distribution  at  about  9=  167  deg.  The 
shock  is  captured,  and  so  it  is  smeared  over  a  mesh  interval. 
T^e  figure  also  shows  the  numerical  results  of  Ref.  IS.  They 
are  also  a  solution  to  Euler's  equations  but  with  the  primary 
cross-flow  shock  was  captured.  In  the  results  of  Ref.  15,  the 
separation  point  was  forced  to  occur  at  9  =  120  deg  in  order  to 
match  the  boundary-layer  separation  point  found  ex¬ 
perimentally.  The  forced  separation  model  of  Ref.  15  is  in 
contradiction  to  the  analytical  work  of  Smith. In  addition, 
there  are  anomalies  in  the  results  of  Ref.  15.  Basically  the 
comparison  of  Fig.  15  shows  that  the  two  results  are  very 
close,  while  the  separation  points  are  very  different.  It  seems 
that  forcing  separation  at  9=  120  deg  simply  inserted  a  wiggle 
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Fig.  15  Surface  pressure  comparison  (<Vf,  =3.  i  =  9.46  deg,  a  =  25 
deg). 
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Fig.  16  Surface  pressure  comparison  ( M „  =  4.25,  3*5  deg.  o  =  25 
deg). 
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in  ihe  surface  pressure  distribution  of  Ref.  15.  The  flow  then 
came  back  to  the  shock-induced  flowfield.  The  cross-flow 
shock  locations  are  very  close.  The  shock  of  Ref.  15  is 
smeared  over  a  few  points  but  is  close  to  the  fit  shock  of  the 
present  calculations.  If  separation  did  occur  at  120  deg,  there 
must  be  a  cross-flow  shock  before  #=  120  deg.  However,  there 
is  no  evidence  of  a  shock  before  separation  in  the  results  of 
Ref.  15. 

Figure  16  shows  the  pressure  distribution  on  a  5-deg  cone  at 
.W*  =4.25  and  a  -  12.35  deg.  Compared  are  the  results  of  the 
present  calculation,  those  computed  using  the  full  potential 
equation  supplied  by  Siclari,  and  the  experimental  results  of 
Rainbird.-'  The  potential  results  of  Siclari  have  been 
corrected  for  nonisentropic  bow  shock  effects  while  main¬ 
taining  the  irroialional  assumption.  Figure  16  indicates  that 
the  Euler  and  potential  results  are  virtually  identical  until  the 
shock.  The  potential  result  docs  not  exhibit  the  minimum  in 
pressure  behind  the  shock  typical  of  separation.  Potential 
calculations  cannot  predict  separation  or  spiraling  without  a 
shedding  sheet.  The  comparison  between  Euler  and  potential 
calculations  affirms  the  fact  that  roiaiionality  is  important 
only  after  the  cross-flow  shock.  A  comparison  between  the 
present  calculation  and  the  experimental  results  of  Rainbird 
clearly  shows  that  the  vorticity  produced  by  the  shock  does 
not  separate  the  flow  near  the  viscous  separation  point.  The 
experimental  data  show  two  separations — primary  at 
(1=120.3  deg  and  secondary  at  160  deg,  while  the  inviscid 
separation  is  at  8=  151.3  deg.  The  longer  plateau  in  pressure 
in  the  experiment  behind  the  primary  separation  point  is  due 
to  the  secondary  separation;  otherwise,  the  Euler  and  ex¬ 
perimental  pressure  distributions  would  be  similar.  It  is  the 
expansion  and  recompression  of  the  reverse  flow  that  causes 
secondary  separation.  It  should  be  clear  that  while  the  vor- 
licity  produced  by  the  cross-tlow  shock  is  not  the  whole 
separation  story,  it  may  play  an  important  role  in  the  process. 
There  are  two  sources  of  vorticity  in  this  flowfield— the  shock 
and  the  boundary  layer.  Both  of  these  sources  of  vorticity 
play  a  role  m  separation  and  the  resulting  spiral  or  vortex. 


Summary  of  Findings  and  Future  Work 

This  study  has  found  that  the  vorticity  produced  by  a  cross- 
flow  shock  present  in  the  supersonic  flow  about  cones  at  high 
incidence  can  cause  separation  and  spiraling.  A  number  of 
important  features  of  these  solutions  were  shown. 

1 1  An  independence  of  these  solutions  to  any  numerical 
artifices;  that  is,  artificial  viscosity  (truncation  error!,  cap¬ 
tured  shocks,  or  initial  conditions. 

2|  Why  the  conceptual  problem  associated  with  the  two- 
dimensional  closed-separation  bubbles  does  not  exist  in 
conical  flow.  Since  no  closed-recirculation  zones  exist  in  the 
conical  flow,  all  quantities  transmitted  along  streamlines  can 
be  uniquely  defined. 

3)  The  variation  of  shock-induced  separation  location  with 
-.hock  strength  (Fig.  12),  indicating  the  possibility  that  the 
spiral  IS  eliminated  before  the  cross-flow  shock. 

4)  The  development  of  a  second  reverse  cross-flow  shock 
(Figs.  14  and  15)  with  the  possibility  of  shock-induced 
secondary  separation. 

It  should  be  obvious  that  any  computation  that  relies  on  the 
potential  How  approximation  alone  cannot  predict  these 
flowfields  While  the  shock-induced  vorticity  may  play  an 
important  role  in  viscous  separation,  it  is  by  no  means  the 
whole  picture  (see  Fig.  I6|. 

This  work  has  uncovered  as  many  problems  as  it  has  solved. 
The  somewhat  surprising  behavior  of  the  lee  plane  sad¬ 
dle  node  combination  once  ihe  spiral  appears  must  be  studied 
in  detail  The  saddle  node  combination  (Figs.  4  and  101  that  is 
present  at  subcritical  »'s  and  remains  with  weak  spirals  is 
eliminated  as  the  spiral  gels  larger  (Figs  4  and  5|  The 


possibility  of  a  second  reverse  shock  (Fig.  15)  causing 
secondary  separation  is  an  interesting  possibility  that  must  be 
investigated. 
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Conical,  Separated  Flows 
with  Shock  and  Shed  Vorticity 

F  Marconi’ 

Grumman  Corporate  Research  Center 
Bethpage,  New  York 

Iniroduciion 

HIS  investigation  deals  with  the  highly  vortical  flow 
about  circular  cones  at  supersonic  speeds  and  high  angle 
of  attack.  The  author  has  shown  that  the  shock  system  of 
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these  Hows  can  cause  separation  and  a  spiral  vortex  near  the 
lee  plane  of  the  cone.'  Earlier,  Salas-  showed  two-dimensional 
Euler  solutions  where  shock  vorticity  had  also  induced  separa¬ 
tion.  The  numerical  investigation  of  Ref.  1  showed  ihat  the 
crossflow  shock  can  produce  large  entropy  gradients  and  vor¬ 
ticity  sufficient  to  cause  separation  and  a  spiral  vortex.  In  ad¬ 
dition  to  shock  vorticitv,  a  separating  boundary  layer  can  shed 
vorticity  into  an  otherwise  irrotational  flow.  This  pheno¬ 
menon  has  been  studied  with  innscid  models  for  many  years. 
Smith'  developed  a  model  for  she  Jding  vorticity  from  the  sur¬ 
face  of  a  smooth  body  into  the  llowfield.  The  present  author'' 
was  ihe  first  lo  use  ihis  model  in  conjunction  with  the  Euler 
equations  to  shed  vorticity  from  both  primary  and  secondary 
crossflow  separation  points.  An  extension  of  the  vorticity 
shedding  model  of  Ref.  3  to  the  situation  of  supersonic 
crossflow  is  used  in  the  present  investigation.  In  the  case  of 
shock-induced  separation,  the  solution  to  the  Euler  equation  is 
unique  in  that  the  separation  point  location  is  computed  along 
with  the  shock  system  and  flowfield.  On  the  other  hand,  the 
separation  point  must  be  presertbed  in  the  case  of  shed  vortic- 
iiy.  There  exists  a  range  of  separation  point  locations  cor¬ 
responding  to  varying  amounts  of  vorticity  being  sheq.  This 
Note  presents  the  results  of  an  investigation  of  the  relationship 
between  shock-vorticity-produced  separation  and  shed- 
voriicity-produced  separation. 

Discussion 

This  study  was  conducted  by  numerically  evaluating  Euler 
solutions  for  the  flow  about  a  5  deg  cone  varying  the  specified 
separation  point  location.  The  freestream  conditions  were 
fixed  at  M,  =4.25  and  oi=  12.35  deg.  The  flow  is  conical  so 
that  only  the  crossflow  plane  (i.e.,  any  plane  normal  to  the 
cone  axis)  need  be  considered.  This  plane  intersects  the  separa¬ 
tion  line  at  a  point  which  will  be  referred  to  as  the  separation 
point.  At  the  large  angle  of  attack  considered  here  the  flow  in 
the  crossflow  plane  is  supercritical  and  so  a  crossflow  shock  is 
present.  With  no  vorticity  shed  from  the  cone  surface,  this 
crossflow  shock  Is  strong  enough  to  produce  enough  vorticity 
to  cause  separation.  In  addition,  shedding  vorticity  from  the 
cone  surface  does  not  eliminate  the  crossflow  shock  so  that 
both  sources  of  vorticity  are  present.  Figures  I  and  2  show  the 
crossflow  streamlines  and  isobars,  respectively,  for  the  flow 
with  no  vorticity  shed  from  the  body,  ^paration  for  this  case 
IS  computed  to  be  at  9,  =  151.3  deg  from  the  wind  plane. 
Figure  1  indicates  that  the  separating  streamlines  leave  the  sur¬ 
face  at  a  large  angle  (57  deg  in  this  particular  case)  relative  to 
the  surface.  When  only  shock  vorticity  is  present  there  is  no 
jump  in  crossflow  velocity  at  the  separation  point,  consistent 
with  the  fact  that  no  vorticity  ;»  being  shed  from  the  surface 
In  the  case  of  the  shock  vorticity  alone,  the  crosvlow 
stagnates  on  both  sides  of  the  separating  streamline,  and  this 
streamline  leaves  the  surface  at  a  large  angle  relative  'o 
Smith'  proved  that  in  order  lo  shed  vorticitv  nto  an  .nherwive 
irrotational  flow  there  must  be  a  jump  in  crossfiow  -eioci  v  .i 
the  separation  point;  the  crossflow  stagnates  -miv  vr  nr  rr 
side  of  the  separating  steamline  It  should  be  pvvinieo  v.  -  i 
in  the  computational  results  presented  here  ai  •  "  - 

shocks  are  captured.  Figure  2  indicates  that  -be  -r-.w  i, 
tured  very  sharply  isee  the  closelv  spaced  v.^ba-.  \.i.; 

ally,  these  captured  shock  results  .omparc  we  ’  ’  r  •  . 

fit  results  of  Ref  1 

The  model  used  to  shed  vo-'i.  -v  "..rr  ■  . 

follows  the  work  o(  Smuh  The  -n.vdr  ■  ,r  -  - 
compuiational  scheme  bv  using  .i  i.-uMr  , 
separation  point  Onegridixuni  ^avcrnr,- 
the  vortex  sheet  and  Ihe  'Ihe'  .  jvN.,'Tir. 
the  case  of  subsoni.  .'.'vvil..^  nr  .; 
of  the  sheet  are  .ompu'rvJ  .i-  r  .  . 
onlv  the  flvvw  'angen.  v  .-.j  -.i"  '  n. 

stagnate  .m  'he  isisos  .jr  .  ,*■.«•.  * 

'tom  'hr  n.'w  vidr  ■  r--  , 
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body  tangentially.  If  the  crossflow  approaching  the  separation 
point  is  supersonic  the  model  used  must  be  changed.  Smith’ 
argued  that  if  the  separating  streamline  leaves  a  smooth  sur¬ 
face  at  a  finite  angle,  then  the  flow  must  stagnate  on  both  sides 
of  the  separation  point  and  no  vonicity  will  be  shed  into  the 
flowfield.  This  argument  is  true  if  the  flow  approaching  the 
separation  point  is  subsonic.  However,  when  the  separating 
flow  is  supersonic  and  the  separating  sheet  leaves  the  body  at 
an  angle  relative  to  it,  then  there  can  be  a  jump  in  crossflow 
velocity  (i.e.,  vorticity  shed).  In  a  locally  supersonic  flow  a 
slope  discontinuity  results  in  a  shock,  not  a  stagnation  point  as 
in  the  subsonic  case.  The  new  model  involves  an  oblique  shock 
at  the  separation  point  unless  the  deflection  is  too  large,  in 
which  case  a  normal  shock  is  formed  and  the  flow  is  separated 
in  a  locally  subsonic  flow.  The  modifications  to  the  computa¬ 
tion  are  as  follows;  the  flow  side  of  the  sheet  is  also  on  the  low- 
pressure  side  of  the  oblique  shock  at  separation  so  that  the 
pressure  there  cannot  be  used  on  the  body  side  of  the  sheet. 
The  crossflow  on  the  body  side  of  the  sheet  is  still  stagnated 
with  the  pressure  there  being  extrapolated  from  the  lee  side  of 
separation.  With  the  conditions  at  the  double  point  defined, 
no  differences  are  taken  across  the  separation  point  on  the 
body.  Once  the  sheet  leaves  the  body  no  provisions  are  made 
for  it  and  it  is  captured  over  a  few  mesh  points.  Reference  4 
gives  more  detail  on  this  procedure. 

Figures  3  and  4  show  the  results  for  the  most  windward 
separation  point  location  computed  (9,  >  115  deg).  Figure  3 
shows  the  crossflow  streamlines.  The  secondary  separation 


Fig.  I  CroMflo*  strtaarihwt  oa  the  S  deg  circaiar  coac  a4.IS, 
a  -  12.35  deg)  scparaiioa  dac  to  slrack  voftletty  aloac  at  f  » 151.3  deg. 


Fig.  2  Isobars  oa  Ike  5  deg  drtalar  coac  (M.  •  4.25,  a  » 12.35  deg) 
scparaiioa  doc  lo  shock  vortidly  alone  al  ISI.3  deg. 


shown  is  due  to  a  strong  reverse  crossflow  shock  (see  Fig.  4); 
note  the  third  vortex  off  the  surface  due  to  the  Interaction  of 
the  two  separations.  In  Fig.  4,  the  isobars  are  shown  and  they 
indicate  an  oblique  shock  at  the  specified  primary  separation 
pent.  The  shock  is  due  to  the  separating  sheet  leaving  the  cone 
surface  at  an  angle  relative  to  it.  A  similar  shock/ sheet  con¬ 
figuration  has  been  noted  experimentally.'  In  the  situation  of 
Figs.  3  and  4.  the  supersonic  crosstlow  passes  through  the 
oblique  shock  and  moves  outward  along  the  separating 
streamline.  The  flow  on  the  lee  side  of  the  sheet  stagnates  at 
the  separation  point.  With  separation  specified  at  5  =  1 15  deg, 
the  jump  in  crossflow  velocity  at  the  separation  point  indicates 
that  significant  vorticity  is  being  shed  from  the  surface.  A 
comparison  of  Figs.  I  and  3  shows  that  the  extent  of  the  vor¬ 
tical  regions  are  comparable,  whereas  the  two  sources  of  vor¬ 
ticity  are  very  different. 

The  relationship  between  shock  vorticity  and  shed  vorticity 
for  conical  separated  flow  is  made  clear  by  considering  Fig.  5. 
The  figure  shows  the  jump  in  crossflow  velocity  (liflcfl 
separation  point  location.  The  jump  in  crossflow  velocity  is 
directly  related  to  the  vorticity  shed  into  the  flowfield  from  the 
separation  point.  Figure  5  represents  all  possible  Euler  solu¬ 
tions  for  the  flow  about  this  cone  at  these  freestream  condi¬ 
tions  and  one  specified  separation  point.  No  solutions  could 
be  found  for  9,  >  151 .3  deg  or  9,  <  1 15  deg.  The  experimental 


Fig.  3  Craraflow  MmuBliacs  oa  the  5  deg  cirealar  coac  (M.  •  4.25, 
12.35  4cg)  segaradoa  forced  al  I  -  115  deg. 


Fig.  4  Isoban  oa  the  5  deg  drcalar  coae  (M.  -  4.25,  a  -  12.35  deg) 
•egaralloa  for  ,  el  f  •  115  deg. 
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data  of  Ref.  S  indicate  primary  separation  at  120  deg  and 
secondary  separation  at  160  deg  for  this  case.  The  shock  con¬ 
figuration  transition  from  an  oblique  crossflow  shock  to  a 
detached  normal  crossflow  shock  occurs  at  about  0, «  128  deg 
(indicated  by  the  shaded  area  in  Fig.  5).  It  should  be  pointed 
out  that  the  jump  in  velocity  at  the  separation  point  in  the 
oblique  shock  cases  was  computed  by  subtracting  the  oblique 
shock  velocity  jump  from  the  numerical  results.  Thus,  the 
jumps  in  velocity  in  Fig.  5  represent  the  jumps  across  the 
vortex  sheet  at  separation.  The  figure  shows  that  this  velocity 
jump  goes  to  zero  smoothly  as  the  separation  point  location 
due  to  shock  vorticity  alone  is  approached  (0,  <=  151.3  deg). 

Conclusion 

The  results  of  this  analysis  indicate  that  there  is  a  relation¬ 
ship  between  separation  produced  by  shock  vorticity  and  shed 
vorticity  and  that  both  sources  of  vorticity  may  be  important. 
The  relationship  is  inferred  because  the  jump  in  velocity  (i.e., 
vorticity  shed)  smoothly  goes  to  zero  as  the  shock  vorticity 
separation  point  location  is  approached.  In  fact,  it  would  seem 
that  separation  due  to  shock  vorticity  alone  can  be  considered 
a  limiting  solution  of  the  set  of  solutions  in  which  vorticity  is 
shed  from  the  surface,  in  this  limiting  solution  the  vonicity 
shed  is  zero. 
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Abstract 

This  work  represents  an  application  of  Prof.  Morettl's  computational 
procedures  to  an  Investigation  of  fluid  physics.  Specifically,  the  effects  of 
shock  vorticity  on  the  supersonic,  inviscid  flow  about  circular,  cross- 
sectional  cones  Is  considered.  While  the  Euler  equations  are  solved,  so  that 
no  viscous  effects  are  Included,  the  shock  vorticity  can  cause  flow  separation 
and  vortex  formation.  In  the  computational  results  to  be  shown  here  the 
Interaction  of  these  vortices  and  shocks  produces  a  number  of  interesting 
phenomena  which  will  be  discussed  In  detail. 

1.  INTROOUaiON 

When  flows  with  significant  vorticity  were  Investigated  by  Salas  [1]  with 
computational  methods,  the  question  of  whether  the  inviscid  flow  separation 
predicted  was  a  numerical  artifice  was  a  serious  concern.  After  the  careful 
shock  fitting  procedures  proposed  by  Moretti  [2]  were  applied  to  an 
Investigation  of  this  phenomenon  [3],  some  of  the  doubts  were  dispelled. 

After  a  number  of  investigations  in  the  area,  it  is  now  accepted  that  shock 
systems  can  produce  enough  vorticity  to  cause  flow  separation.  While  Salas 
[1]  may  have  been  the  first  researcher  In  modern  computational  fluid  dynamics 
to  investigate  this  phenomenon,  Kuchemann  [4]  considered  the  possibility  of 
shock  vorticity  causing  separation  In  1967.  In  addition,  Frankel  [5] 
developed  analytical  solutions  to  the  Euler  equations  which  exhibited 
separation  due  to  distributed  vorticity.  The  current  paper  will  show  that  the 
shock  system  produced  by  a  circular  cone  flying  at  supersonic  speeds  and  high 
angle  of  attack  can  generate  enough  distributed  vorticity  to  cause  separation 
and  vortex  formation. 
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Figure  1  shows  a  sketch  of  the  flow  field  that  is  considered  here.  The 
geometry  is  a  slender  cone  that  is  immersed  in  a  high  speed  flow  at  a  high 
angle  of  attack.  As  the  angle  of  attack  is  increased  past  the  cone  half 
angle,  the  crossflow  component  tangent  to  the  body  (u)  will  become 
supersonic.  Due  to  symmetry,  this  velocity  must  go  to  zero  in  the  lee  plane, 
causing  a  crossflow  shock  to  develop  on  the  lee  side  of  the  body.  As  the 
angle  of  attack  is  raised,  this  shock's  strength  increases.  The  strongest 
portion  of  this  shock  is  at  the  cone  surface,  and  it  diminishes  as  it  moves 
into  the  field.  It  goes  to  zero  strength  relatively  close  to  the  body.  It  is 
this  rapid  decrease  in  shock  strength  which  produces  enough  vorticity  to  cause 
vortex  formation.  It  was  Smith  [6]  who  showed  that  in  order  to  have  a  vortex 
in  conical  flow  a  non-zero  component  of  radial  vorticity  is  required. 

Figure  1  shows  the  simplest  shock/vortex  pattern  encountered  during  this 
work.  The  interaction  of  the  shock  and  vortex  can  produce  complex  shock 
patterns  and  multiple  vortices.  In  addition,  aside  from  shock  vorticity, 
vorticity  can  be  shed  from  the  surface  of  the  cone.  The  present  author  [7) 
applied  a  model  to  the  Euler  equations  for  shedding  vorticity  from  the  surface 
of  smooth  bodies.  In  Ref.  [8]  the  relationship  between  shock  and  shed 
vorticity  is  studied.  In  this  paper  only  shock  vorticity  will  be 
considered.  The  flow  about  a  10*  half  angle  cone  at  two  Mach  numbers,  one 
intermediate  5  and  the  other  low  hypersonic  10,  will  be  analyzed.  As  the 
angle  of  attack  is  increased  from  10*  to  25*,  shock/vortex  patterns  of 
increasing  complexity  result.  Before  these  interesting  flow  fields  are 
discussed,  the  computational  procedures  used  will  be  outlined. 

2.  COMPUTATIONAL  PROCEDURE 

The  fully  three-dimensional  Euler's  equations  are  solved  with  an  explicit 
marching  technique.  The  marching  direction,  z  (Fig.  1),  is  an  iterative 
coordinate  for  the  conical  flow  considered  here.  The  scheme  is  restricted  by 
the  fact  that  the  axial  component  of  velocity,  W  (Fig.  1),  must  be  supersonic 
everywhere.  The  marching  scheme  is  continued  until  the  flow  field  is 
invariant  with  respect  to  the  computational  marching  direction  except  for  a 
scale  factor.  The  finite  difference  scheme  used  is  Moretti's  characteristic 
based  x-scheme  [9].  The  bow  shock  and  the  primary  crossflow  shock  are  fit  and 
are  forced  to  satisfy  the  exact  Rankine-Hugoniot  Jump  conditions.  The  bow 
shock  is  fit  as  the  outermost  boundary  of  the  flow  field.  On  the  low  pressure 


side  of  the  bow  shock,  freestream  conditions  exist.  The  crossflow  shock  Is 
fit  as  an  Internal  boundary  of  the  flow  field  with  Its  low  pressure  side  being 
computed  as  the  computation  proceeds  [10]. 

The  crossflow  shock  computation  Is  a  critical  part  of  the  overall 
procedure,  particularly  since  It  plays  such  a  critical  role  In  separation.  As 
mentioned  In  the  Introduction,  It  Is  the  vorticity  produced  by  the  crossflow 
shock  which  Is  significant.  Any  scheme  that  captures  the  crossflow  shock  and 
Introduces  additional  artificial  viscosity  to  stabilize  It  runs  the  risk  of 
distorting  the  separation  and  vortex.  The  bow  shock  also  Introduces  vorticity 
Into  the  flow  field  but  not  enough  to  produce  separation.  The  low  pressure 
side  of  the  crossflow  shock  Is  computed  with  one-sided  differences  away  from 
the  shock.  This  Is  consistent  with  the  fact  that  the  crossflow  Is 
supersonic.  The  Rankine-Hugonlot  conditions,  together  with  the  compatibility 
condition  along  a  bicharacteristic  reaching  the  shock  on  Its  high  pressure 
side,  supply  enough  Information  to  compute  the  deviation  of  each  shock  point 
from  a  conical  ray  In  addition  to  all  the  primitive  variables  on  Its  high 
pressure  side.  The  bicharacteristic  used  Is  the  one  In  the  plane  containing 
the  local  normal  to  the  shock  and  the  marching  direction.  All  shock  points 
are  computed  with  the  post  correction  scheme  proposed  by  Rudman  [11]  and 
Independently  by  deNeef  [12] .  The  shocks  converge  (I.e.,  each  shock  point 
becomes  aligned  with  a  conical  ray)  with  the  rest  of  the  flow  field.  In  all 
the  computations  presented  here,  the  last  crossflow  shock  point  fit  had  a 
normal  Mach  number  of  approximately  1.05  (pressure  ratio  of  1.12).  The  finite 
difference  scheme  was  able  to  capture  weaker  shock  points.  In  addition,  the 
finite  difference  scheme  was  used  to  capture  all  reverse  crossflow  shocks  and 
oblique  crossflow  shocks  In  the  flow  field  (to  be  discussed  later).  These 
shocks  are  usually  weak  and  the  scheme  can  capture  them  accurately. 

The  computational  grid  used  here  was  developed  from  the  polar  coordinate 
system  (r,e,z)  shown  In  Fig.  1.  The  grid  Is  clustered  exponentially  In  the 
radial  (r)  direction  with  a  stretching  of  the  form 

r  •  -in(l  -  (l-e'®)X]  (C-B)/6+B  (1) 

where  X  Is  the  computational  coordinate  going  from  0  on  the  body  to  1  at  the 
bow  shock.  The  other  computational  coordinate  Is  Y,  which  Is  a  linearly 
normalized  polar  angle  e.  Y  Is  normalized  to  go  from  0  In  the  windward 
symmetry  plane  to  1  on  the  low  pressure  side  of  the  crossflow  shock  and  from  0 


again  on  the  high  pressure  side  of  the  shock  to  1  In  the  lee  plane.  The 
transformation  to  the  computational  plane  Is  completed  by  Z  »  2.  In  the 
transformation,  Eg  (1),  B  and  C  are  the  polar  radii  of  the  body  and  bow  shock, 
and  B  Is  a  clustering  parameter.  The  grid  points  are  evenly  spaced  In  each 
X-Y  plane.  The  transformation  clusters  grid  points  near  the  body  depending  on 
the  value  of  s  used,  and  It  approaches  linear  as  s  goes  to  zero.  All  the 
calculations  shown  here  were  performed  with  a  *  2.5.  The  gradient  In  the 
strength  of  the  crossflow  shock  Is  quite  large,  which  Is  the  reason  for  the 
large  vorticity  produced  behind  It.  The  clustering  of  grid  points  near  the 
body  allows  for  enough  resolution  to  compute  this  shock  decay  properly.  The 
number  of  points  required  to  resolve  this  phenomenon  without  a  stretching 
would  be  Inordinately  large.  In  all  the  calculations  shown  here  the  cross- 
sectional  grid  had  89x89  points.  There  were  anywhere  from  13  to  20  points  on 
the  lee  side  of  the  crossflow  shock  depending  on  Its  position;  the  rest  of  the 
89  points  were  evenly  spaced  on  the  low  pressure  side  of  the  shock. 

Shock  Induced  Inviscid  separation  produces  a  contact  surface  emanating 
from  the  separation  point  (Fig.  1).  It  Is  this  contact  that  ultimately 
spirals  to  form  a  vortex.  The  contact  surface  has  a  jump  In  entropy  and  an 
accompanying  jump  1n  velocity.  In  conical  flow,  entropy  Is  constant  on 
crossflow  streamlines.  Thus,  one  can  see  (F1g.  1)  that  the  entropy  on  the 
windward  side  of  the  separating  contact  comes  from  the  windward  stagnation 
point  of  the  cone.  The  streamline  that  wets  the  body  passes  through  the  base 
of  the  crossflow  shock  to  form  the  high  entropy  side  of  the  contact.  There 
are  a  number  of  possibilities  for  the  entropy  on  the  lee  side  of  the 
contact.  The  crossflow  streamlines  that  wet  the  body  are  tracked  In  each  step 
of  the  Iteration,  and  the  proper  entropy  Is  Imposed  on  the  surface  including 
the  jump  In  entropy  at  the  separation  point.  All  entropy  discontinuity  off 
the  body  Is  captured  with  the  finite  difference  scheme. 

The  finite  difference  scheme  used  1n  this  work  Is  an  explicit  marching 
scheme,  which  Is  notoriously  Inefficient  for  converging  to  a  conical  or  steady 
solution.  One  advantage  of  the  explicit  marching  scheme  Is  that  It  Is  totally 
"vectorizable."  The  time  consuming  parts  of  the  code  (I.e.,  Interior  point 
computation)  utilized  the  vector  architecture  of  the  Cray  XMP  computer.  The 
computations  shown  here  typically  took  45  CPU  minutes  on  the  Cray  XMP.  In 


each  case  shown,  the  maximum  residual  (I.e.,  the  derivative  of  pressure  In  the 
marching  direction)  was  reduced  at  least  six  orders  of  magnitude. 

3.  COMPLEX  SHOCK/VORTEX  PATTERNS 

In  this  computational  investigation,  a  single  cone  Is  considered.  The 
half  angle  1s  10”,  and  two  Mach  numbers  (5  and  10)  were  studied.  At  both  Mach 
numbers  the  angle  of  attack  was  Increased  from  10"  to  25"  In  Increments  of 
5".  Of  Interest  Is  the  development  of  the  shock/vortex  pattern  with 
Increasing  angle  of  attack. 

First  the  low  Mach  number  (5)  will  be  discussed.  When  the  angle  of 
attack  Is  equal  to  the  cone  angle  of  10",  the  cross  flow  does  not  expand 
enough  to  become  supercritical.  The  maximum  crossflow  Mach  number  on  the 
surface  Is  only  0.93.  Figure  2  shows  the  crossflow  streamline  pattern  and  the 
shock  system  (note  the  absent  of  a  crossflow  shock).  The  crossflow 
streamlines  are  traced  by  projecting  the  conical  velocity  components  (u  and  v 
of  Fig.  1)  onto  a  plane.  If  the  angle  of  attack  Is  Increased  by  5"  to  15°,  a 
small  region  of  supersonic  crossflow  does  appear  In  the  flow  field.  Figure  3 
shows  the  sonic  bubble  and  an  accompanying  crossflow  shock.  Figure  4  shows 
the  crossflow  streamlines  for  this  case.  Figure  4a  shows  the  whole  cross- 
sectional  flow  field  while  Fig.  4b  shows  only  the  region  near  the  lee  plane. 

In  Fig.  4b  a  small  vortex  Is  shown  very  near  the  lee  plane.  The  crossflow 
shock  has  a  maximum  normal  Mach  number  of  1.66  at  the  body.  This  shock  is 
strong  enough  to  produce  only  a  very  small  vortex. 

When  the  angle  of  attack  Is  Increased  to  20”,  the  shock  gets 
significantly  stronger  (maximum  shock  Mach  number  of  2.46)  and  the  resulting 
vortex  Is  large.  Figure  5  shows  the  crossflow  streamlines  and  shock  system  at 
a  »  20".  While  the  shock/vortex  system  Is  significant.  It  Is  simple  In  the 
sense  that  there  Is  still  one  vortex  and  one  crossflow  shock.  Figure  6  shows 
the  Isobars  for  this  case,  and  Fig.  7  shows  the  surface  pressure  distribution. 
The  closed  Isobars  of  Fig.  6  represent  a  minimum  In  pressure  at  the  center  of 
the  vortex.  The  surface  pressure  shows  a  local  minimum  (Fig.  7,  e  ■  170") 
directly  under  the  vortex.  The  reverse  crossflow  expands  from  the  lee  plane 
to  achieve  It  maximum  velocity  at  the  minimum  In  pressure  at  a  ■  170*.  The 
flow  then  recompresses  to  go  to  zero  crossflow  velocity  at  the  separation 
point  a  •  165*.  This  recompression  can  become  a  shock  as  will  be  seen  at 
higher  angle  of  attack  when  the  shock  pattern  becomes  complex. 
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Figure  8  shows  contour  lines  of  the  radial  component  of  vorticity: 

Op  »  -  (u  cos  a  +  1^  sin  a  -  sin  a) 

where  u  and  v  are  the  spherical  crossflow  velocities  and  R,  e,  a  are  the 
spherical  coordinates  shown  In  Fig.  1.  Actually,  Fig.  8  shows  Rfl[^  contours, 
which  Is  the  appropriate  conical  quantity.  The  figure  shows  that  the 
vorticity  Is  distributed  somewhat  uniformly  along  the  lower  portion  of  the 
shock.  The  uniform  vorticity  Is  carried  downstream  and  swept  Into  the 
vortex.  This  vorticity  distribution  can  be  understood  more  clearly  by 
considering  the  entropy  distribution  along  the  shock.  Figure  9  shows  the 
entropy  distributions  on  the  low  and  high  pressure  sides  of  the  shock. 
According  to  Crocco  theorem,  the  vorticity  Is  related  to  the  gradient  of 
entropy.  The  highest  gradient  of  entropy  Is  at  the  surface  on  the  high 
pressure  side  of  the  shock.  The  entropy  gradient  decreases  away  from  the 
surface  and  seems  to  have  a  constant  value  for  s/r  >  0.3.  This  (s/r  «  0.3)  Is 
where  the  shock  has  a  rapid  change  In  slope.  Figure  8  Indicates  that  the 
vorticity  Is  constant  after  this  point.  The  rapid  change  In  shock  slope  seems 
to  occur  at  a  radius  equal  to  that  at  the  top  of  the  vortex.  At  even  larger 
angles  of  attack,  this  slope  change  becomes  a  discontinuity  and  the  shock 
system  becomes  complex. 

At  a  freestream  Mach  number  of  5,  the  shock  system  becomes  complex  at  an 
angle  of  attack  of  25*.  Figure  10  shows  the  crossflow  streamlines  and  Fig.  11 
the  Isobars.  The  rapid  change  In  shock  slope  has  become  a  discontinuity  as 
the  shock  extends  over  the  top  of  the  vortex.  This  discontinuity  In  shock 
slope  Is  accompanied  by  an  additional  shock  which  extends  toward  the  lee 
plane.  This  shock  Is  captured  but  Is  clearly  seen  In  the  Isobar  plot  of  F1g. 
11.  The  crossflow  shock  Is  no  longer  simple  but  has  a  triple  point  separating 
two  portions  of  shock. 

Around  the  primary  vortex  a  "pin  wheel"  (13)  shock  pattern  Is  formed. 

The  crossflow  shock  becomes  oblique  to  the  flow  Just  off  the  surface  so  that 
the  cross  flow  behind  the  shock  Is  supersonic  Just  off  the  surface;  see  the 
sonic  line  plot  of  Fig.  12.  The  bundle  of  supersonic  cross  flow  extends  to 
the  triple  point.  This  flow  must  be  shocked  In  order  to  meet  the  symmetry 
condition  In  the  lee  plane,  and  so  a  crossflow  shock  sits  on  top  of  the  vortex 
(Fig.  11).  As  lee  plane  flow  Is  accelerated  by  the  vortex.  It  becomes 
supersonic  so  that  another  sonic  bubble  Is  formed  attached  to  the  lee  plane 


(Fig.  12).  This  supersonic  flow  must  be  shocked  to  meet  the  surface  boundary 
condition.  Note  In  Fig.  11  that  the  shock  Is  normal  to  the  lee  plane  just 
above  the  body.  Finally,  the  surface  flow  Is  accelerated  under  the  primary 
vortex  and  becomes  supersonic  (see  the  third  sonic  bubble  In  Fig.  12).  This 
flow  Is  shocked  In  order  to  become  tangent  to  the  Inside  of  the  separation 
line.  This  process  of  expanding  and  shocking  as  the  flow  moves  around  the 
vortex  sets  up  a  system  of  three  shocks  which  the  author  first  referred  to  as 
the  pin  wheel  type  In  Ref.  13.  This  shock  system  Is  clearly  shown  In  the 
Isobars  of  Fig.  11. 

Figure  10  shows  a  second  small  vortex  just  off  the  lee  plane  above  the 
primary  vortex.  This  Is  a  spiral  node  caused.  It  seems,  by  the  vortlcity 
produced  above  the  triple  point  on  the  main  crossflow  shock.  A  plot  of 
entropy  contours  (Fig.  13)  Indicates  that  the  contact  surfaces  at  the  separa¬ 
tion  point  and  the  triple  point  are  captured  quite  sharply  by  the  finite 
difference  scheme.  Figure  13  also  Indicates  that  the  captured  shocks  have  no 
entropy  Jumps  associated  with  them.  This  Inaccuracy  of  the  scheme  Is  minimal 
since  the  normal  Mach  number  to  any  of  these  shocks  Is  small  (I.e.,  maximum  of 
1.52).  Figure  13  also  shows  the  entropy  singularity  on  the  lee  plane  just 
above  the  secondary  vortex.  This  Is  due  to  the  coalescence  of  streamlines 
with  vastly  different  entropy.  For  the  sake  of  comparison  the  surface  pres¬ 
sure  distribution  Is  shown  In  Fig.  14.  This  figure  should  be  compared  to  Fig. 
7  In  which  the  angle  of  attack  Is  20*.  The  comparison  Indicates  how  much 
stronger  the  crossflow  shock  has  become  with  a  5*  Increase  In  angle  of  attack 
and  how  the  recompression  under  the  vortex  has  turned  Into  a  reverse  crossflow 
shock. 

Ihe  same  10*  half  angle  cone  was  Investigated  In  the  low  hypersonic 
regime.  At  Mach  10  the  10*  angle  of  attack  case  was  attached  as  was  the 
M,  ■  5,  o  ■  10*  case.  Even  though  at  Mach  10  the  crossflow  Is  supersonic,  no 
crossflow  shock  was  present.  Figure  15  compares  the  crossflow  Mach  numbers  at 
M^  ■  5  and  10  at  a  >  10*.  The  figure  shows  that  the  M.  ■  10  case  has  a 
slightly  supercritical  bubble  (maximum  level  Mach  number  of  1.15)  but  not  high 
enough  to  form  a  shock  or  a  vortex.  As  the  angle  of  attack  Is  Increased  by  5° 
to  15*,  a  crossflow  shock  and  vortex  do  form.  Figure  16  shows  the  crossflow 
streamlines  and  the  shock.  In  this  case  the  shock  Is  significant  (maximum 
normal  Mach  number  of  2),  but  the  resulting  vortex  Is  small. 
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Complex  shock  patterns  occur  at  a  =  20"  for  =  10  Instead  of  a  =  25°  at 
»  5.  Figures  17  and  18  show  the  crossflow  streamlines  and  Isobars, 
respectively,  for  the  *  10  and  a  *  20“  case.  Figure  18  shows  the  pin  wheel 
type  shock  pattern  discussed  previously  In  addition  to  the  triple  point  in  the 
crossflow  shock.  It  should  be  noted  that  there  Is  no  secondary  vortex  in  this 
case;  Instead,  there  Is  a  simple  node  in  the  lee  plane  above  the  primary 
vortex. 

The  last  case  considered  Is  »  10  and  a  *  25°.  Figure  19  shows  the 
streamlines  and  Fig.  20  the  Isobars.  Again  the  pin  wheel  shock  pattern  is 
quite  evident  In  addition  to  the  triple  point.  The  primary  crossflow  shock 
has  a  maximum  normal  Hach  number  of  3.8  at  the  surface  and  Is  quite  strong. 

The  captured  reverse  crossflow  shock  has  a  normal  Mach  number  of  1.8.  This 
has  a  secondary  vortex  above  the  main  vortex  In  addition  to  a  simple  node  In 
the  lee  plane  halfway  between  the  body  and  the  bow  shock  (not  shown  In  the 
figure).  The  secondary  spiral  node  above  the  main  vortex  Is  seen  more  clearly 
In  the  expanded  view  of  Its  region  shown  In  Fig.  21.  Figure  21  Is  an  expanded 
view  of  the  dashed  box  In  Fig.  19.  It  can  be  seen  that  the  node  does  not  form 
a  large  vortex  as  does  the  primary  vortex  but  It  does  seem  to  spiral  at  the 
node. 
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FIGURE  LEGENDS 


1.  Sketch  of  flow  field  and  coordinate  systems. 

2.  Crossflow  streamlines  and  shocks  10*  cone,  M  =  5,  a  =  10°. 

3.  Crossflow  sonic  line  and  shock  10°  cone,  =  5,  a  =  15°. 

4.  Crossflow  streamlines  and  shocks. 

a.  entire  cross  section 

b.  region  near  lee  plane 

5.  Crossflow  streamlines  and  shock  10°  cone,  *  5,  o  «  20°. 

6.  Cross  sectional  Isobars  and  shock  10°  cone,  =  5,  o  =  20°. 

7.  Surface  pressure  distribution  10°  cone,  =  5,  a  =  20°. 

8.  Cross  sectional  vorticity  contours  10°  cone,  «  5,  a  =  20° 

9.  Entropy  distribution  along  the  crossflow  shock  10°  cone, 

(s  In  the  arch  length  above  the  shock,  r  Is  the  cone  radius) 

10.  Crossflow  streamlines  and  shock  10°  cone,  «  5,  o  »  25°. 

11.  Isobars  and  shocks  10°  cone,  «  5,  o  »  25°. 

12.  Crossflow  sonic  lines  and  shock  10°  cone,  *  5,  o  »  25°. 

13.  Cross  sectional  entropy  contours  and  shock  10°  cone,  =  5, 

14.  Surface  pressure  distribution  10°  cone,  =  5,  o  =  25°. 

15.  Surface  Mach  number  comparison  10°  cone,  o  »  10°,  *  5  and 

16.  Crossflow  streamlines  and  shocks  10°  cone,  M_  ■  10,  a  *  15°. 

17.  Crossflow  streamlines  and  shock  10°  cone,  »  10,  o  »  20°. 

18.  Isobars  and  shocks  10°  cone,  M^  *  10,  o  »  20°. 

19.  Crossflow  streamlines  and  shocks  10°  cone,  M  »  10,  o  *  25°. 

20.  Isobars  and  shocks  10°  cone,  M^  »  10,  a  «  25°. 

21.  Expanded  view  of  crossflow  streamline  10°  cone,  »  10,  a  ■ 


